NPS71-87-001   (REVISED) 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


NOTES  ON  SEARCH,  DETECTION 
AND  LOCALIZATION  MODELING 


R.  N.  FORREST 


APRIL  1987 
(Revised  October  1987) 


Approved  for  public  release;  distribution  unlimited, 

Prepared  for: 

Naval   Postgraduate  School 

Monterey,  CA  93943-5000 


FedDocs 

D    208.14/2 

NPS-7 1-87-00 1/987 


tedU 
UK 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


Rear  Admiral  R.  C.  Austin  K.  T.  Marshall 

Superintendent  Acting  Provost 

Reproduction  of  all  or  part  of  this  report  is  authorized. 


UNCLASSIFIED 


ISECURITY  CLASSIFICATION  OF  THIS   PAGE 


DUDLEY  KNOX  LIBRARY 

y^  P0ST?RAD_L<ATE  SCHOOL 


REPORT  DOCUMENTATION  PAGE       ""^""^  CA  93943-5101 


1a    REPORT  SECURITY  CLASSIFICATION 


UNCLASSIFIED 


lb    RESTRICTIVE   MARKINGS 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b    DECLASSIFICATION /DOWNGRADING  SCHEDULE 


3     DISTRIBUTION /AVAILABILITY  OF   REPORT 

Approved  for  public  release;  distribution 
unlimited. 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

NPS71-87-001    (REVISED) 


5    MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a    NAME  OF  PERFORMING  ORGANIZATION 

Naval   Postgraduate  School 


6b    OFFICE  SYMBOL 
(If  applicable) 

Code  71 


7a    NAME  OF  MONITORING  ORGANIZATION 


6c.  ADDRESS  (City.  State,  and  ZIP  Code) 

Monterey,  CA  93943-5000 


7b    ADDRESS  (City.  State,  and  ZIP  Code) 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Naval  Postgraduate  School 


8b    OFFICE  SYMBOL 
(If  applicable) 


9    PROCUREMENT  INSTRUMENT  IDENTIFICATION   NUMBER 


8c.  ADDRESS  (City,  State,  and  ZIP  Code) 

Monterey,  CA  93943-5000 


10    SOURCE  OF  FUNDING   NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 
NO 


TASK 
NO 


WORK   UNIT 
ACCESSION   NO 


11.  TITLE  (Include  Security  Classification) 

NOTES  ON  SEARCH,  DETECTION  AND  LOCALIZATION  MODELING 


12    PERSONAL  AUTHOR(S) 


Forrest,   R.   N, 


13a    TYPE  OF  REPORT 

Technical 


13b    TIME  COVERED 
FROM  TO 


14    DATE  OF  REPORT   ( Year,  Month,  Day) 

1987  Oct 


15    PAGE  COUNT 

75 


16   SUPPLEMENTARY  NOTATION 


17 


COSATI  CODES 


FIELD 


GROUP 


SUB-GROUP 


18   SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

search  modeling,  detection  modeling,  localization  modeling 


19   ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

The  report  is  an  extended  outline  of  some  of  the  course  material  used  in 
courses  on  search,  detection  and  localization  modeling. 


20    DISTRIBUTION/AVAILABILITY  OF  ABSTRACT 

CJO  UNCLASSIFIED/UNLIMITED       □  SAME  AS  RPT  Q  DTIC  USERS 


21    ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


22a    NAME  OF  RESPONSIBLE  INDIVIDUAL 

R.  N.   Forrest 


22b  TELEPHONE  (Include  Area  Code) 


(408)646-2653 


22c    OFFICE  SYMBOL 

Code  55Fo 


DDFORM  1473,  84  mar 


83  APR  edition  may  be  used  until  exhausted 
All  other  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

ft  U.S.  Government  Priming  Office:    1986— 60S  243 


Preface 


This  report  is  a  collection  of  material  that  has  been  used 
in  courses  on  search,  detection  and  localization  modeling.   Its 
organization  follows  to  some  extent  material  by  S.  M.  Pollock  in 
Selected  Methods  and  Models  in  Military  Operations  Research  which 
is  listed  in  the  report  bibliography.   The  report  is  not  intended 
to  be  a  text  on  these  subjects.   In  particular,  in  some  areas  it 
does  not  provide  the  depth  of  coverage  that  is  found  in  the  book 
Search  and  Detection  by  Alan  R.  Washburn  which  is  cited  in  this 
report  as  Reference  21. 

Typographical  and  other  errors  that  have  been  found  in  the 
original  version  of  this  report  have  been  corrected  and  several 
sections  have  been  modified  or  expanded  in  this  revision. 
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I.   Detection  Models  and  Signal  Detection  Theory 

Signal  detection  theory  is  the  basis  for  analyzing  the 
detection  models  that  are  described  in  this  report.   In  signal 
detection  theory,  the  decision  making  portion  of  a  detection 
system  is  called  the  receiver  and  a  detection  experiment  is  the 
observation  by  a  receiver  of  input  data  generated  during  some 
time  interval.   The  data  that  is  related  to  a  target  is  called 
signal.   The  data  that  is  not  related  to  the  target  is  called 
noise.   In  general,  the  target  data  is  associated  with  a 
localization  region  that  in  some  cases  is  called  a  resolution 
cell.   When  a  detection  experiment  is  performed,  either  the  event 
H0  =  {the  receiver  input  is  noise}   or  its  complement 
H^  =  {the  receiver  input  is  signal  and  noise)   will  occur.   In 
the  first  detection  models  that  are  described  here,  after 
analysis  of  the  input  data  by  a  receiver,  either  the  event 
Dq  ={the  receiver  decides  the  input  is  noise}   or  its  complement 
D^  =  {the  receiver  decides  the  input  is  signal  and  noise}   will 
occur.   Detection  models  for  which  D-±     is  the  complement  of  Dq 
are  called  binary  detection  models  or  forced  choice  detection 
models.   Four  events  which  are  important  in  binary  detection 
model  are  indicated  in  the  Venn  diagram  of  Figure  1. 

The  Venn  diagram  emphasizes  a  decision  problem  that  is 
associated  with  a  receiver  that  can  be  modeled  using  a  binary 
detection  model.   The  problem  is  this:   Under  what  conditions 
should  the  event  D-±     occur?  That  is,  under  what  conditions 
should  a  receiver  decide  that  the  input  data  indicates  a  target 
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was  present  in  an  observed  localization  region  either  at  the  time 
or  prior  to  the  time  of  the  observation? 


H, 


H- 


D0   n  H0 

D0   n   Hx 

Di   n  H0 

Dx    n   Hx 

Figure  1.   Four  events  of  importance  in  binary  detection  models. 


In  the  detection  model  descriptions  that  follow,  the 
following  notation  and  terminology  are  used:   Pf  =  PfDjjHg),   the 
probability  of   D^  given  H0,  is  called  the  false  alarm 
probability;   p<j  =  P(D1|H1)/  the  probability  of  D-±      given  H^, 
is  called  the  detection  probability  and   P  =  P(Hj_) ,   the 
probability  of  H^,  is  called  the  prior  probability.   It  is  the 
probability  that  a  target  will  be  in  the  localization  region  when 
the  input  data  is  generated. 

In  the  detection  models,  the  input  to  a  receiver  consists  of 
a  sequence  Yii •  *  *Ym     tna^  is  a  realization  of  a  discrete 
parameter  stochastic  process  whose  random  variables  are  elements 
of  a  continuous  parameter  stochastic  process.   Although  the 
receiver  input  in  some  cases  can  be  considered  to  be  described  by 
a  continuous  parameter  stochastic  process,  because  of  the  finite 
amount  of  information  (unique  data)  contained  in  a  bounded 
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sequence  of  finite  length,  a  discrete  parameter  stochastic 
process  is  sufficient  to  represent  the  receiver  input  in  these 
cases.   This  is  established  formally  by  the  stochastic  sampling 
theorem. 

In  the  models,  for  an  observation,  the  input  stochastic 
process  has  the  following  characteristics:   It  is  a  noise  process 
when  there  is  no  target  data  and  it  is  a  noise  process  plus  a 
signal  process  when  there  is  target  data.   In  the  first  detection 
model  that  is  described  in  Section  III,  the  signal  process  is  a 
deterministic  process.   In  the  second  and  third  detection  models 
that  are  described  there,  the  signal  process  is  a  random  process. 
To  specify  the  noise  process  or  a  random  signal  process,  one 
needs  only  to  specify  the  joint  distribution  of  the  finite 
sequence  of  random  variables  that  determine  the  process.   If  the 
signal  process  is  a  deterministic  process,  the  signal  values  are 
a  sequence  of  values  that  are  assumed  to  be  known  before  an 
observation  is  performed.   Consequently,  in  this  case,  to  specify 
the  process,  one  needs  only  to  specify  these  values. 


II.   Decision  Criteria 

To  simplify  the  discussion  of  decision  criteria  and 
decision  rules,  a  receiver's  input  will  be  assumed  to  be 
determined  by  a  single  decision  random  variable  Y.   In  this 
case,  the  input  process  in  determined  by  the  conditional 
distribution  function  Fy(y|H0)   when  the  input  is  noise  alone 
and  by  the  conditional  distribution  function  FyfylH!)   when  the 
input  is  signal  plus  noise. 

The  condition  that  a  receiver's  input  is  required  to 
satisfy  in  order  that  the  event  D-±     will  occur  can  be  specified 
in  terms  of  a  decision  rule.   For  the  assumed  case,  a  decision 
rule  is  a  rule  which  determines  for  every  observable  value  of  Y 
the  decision  that  the  receiver  is  to  make.   The  decision  rule  can 
be  considered  to  be  a  function   0(y)   which  relates  each 
observable  value  y  to  one  or  the  other  of  the  following  two 
commands:   d0  =  "Decide  that  the  receiver  input  was  noise."  and 
d-±   =  "Decide  the  input  was  signal  and  noise.".   Choosing  a 
decision  rule  <p(y)      defines  a  set  n  of  observable  values  of  Y 
such  that  the  event  D±   =  {  Y  c  n  } . 

The  problem  which  was  considered  in  Section  I  can  now  be 
restated  in  the  following  way:   What  criterion  should  be  adopted 
in  order  to  determine  a  decision  rule  or,  equivalently,  its 
corresponding  set  n  ?  A  desirable  characteristic  for  a  criterion 
is  suggested  by  the  following  argument:   Consider  the  odds  in 
favor  of  Hi  given  y   is  observed.   That  is,  consider 


PfHjjY  =  y)/P(H0|Y  =  y)  .   One  might  expect  that  y  would  be  a 
member  of  the  set  n   if  and  only  if  y  made  this  ratio  equal  to 
or  greater  than  some  value  k.   But  this  is  equivalent  to  defining 
n   as  follows:   n  =  {  y  :  L(y)  >  K  }  where  L(y)   is  the 
likelihood  ratio  associated  with  an  observed  value  y  and  K  is 
a  constant  related  to  the  constant  k.   This  suggests  that 
choosing  an  optimum  criterion  is  equivalent  to  choosing  an 
optimum  value  for  K.   Four  specific  decision  criteria  are  defined 
next  in  terms  of  K.   For  each  criterion,   n   has  the  above 
form.   But  for  each  criterion  the  choice  of  K  is  different. 
The  decision  criteria  are: 

1.  The  Neyman-Pearson  Criterion:   Choose  n   so  that  p^ 

is  a  maximum  subject  to  the  constraint  that  pf  <  a  where  a   is 
a  specified  value.   For  a  continuous  decision  random  variable, 
the  constant  K  is  chosen  so  that  pf  =  a. 

2.  The  Bayes  Criterion:   Choose  n   so  that  the  expected 
cost  of  a  receiver's  decision  is  a  minimum.   For  a  continuous 
decision  random  variable,   K  =  [  (c^q-c0o)/  (coi""cll)  1  *  (1_p)/p   if 
c10  >  c00   an<*  c01  >  cll  where  c^j   is  the  cost  of   D^  n  Hj . 

3.  The  Ideal  Observer  Criterion:  Choose  n  so  that  the 
probability  that  the  receiver  makes  an  incorrect  decision  is  a 
minimum.   K  =  (1-P)/P  for  a  continuous  decision  random  variable. 

4.  The  Minimax  Criterion:   Choose  n  when  P  is  unknown  so 
that  the  maximum  expected  cost  of  a  receiver's  decision  is  a 
minimum.   If   c^q  >  c00   anc*   c01  >  cll'  tnen 


K  =  [ (cio~cOo)/ (c01~cll) ] * (1-P*)/P*   for  a  continuous  decision 
random  variable.   Here,   P*  is  the  value  of  the  prior  probability 
P  that  would  make  the  expected  cost  of  a  receiver's  decision  a 
maximum  if  P  were  known  and  the  Bayes  Criterion  were  used. 

If  a  model  which  specifies  the  conditional  distributions 
Fy(ylHo)   and  Fy(ylHl)   an<*  a  decision  rule  are  adopted,  then 
the  value  of  pf  and  the  value  of  p^  are  determined.   This 
pair  of  values   (PfrPd)   ^s  called  a  receiver  operating  point. 
If  the  decision  rule  results  from  using  a  likelihood  ratio 
criterion  such  as  one  of  the  four  listed  above,  then  it  will 
involve  the  parameter  K  since  n  =  {y:  L(y)  >  K) .   And,  for  a 
given  value  of  K,  since  fi  uniquely  determines  the  pair 
(Pf ,pa) ,  a  single  operating  point  results.   By  varying  K,  a  set 
of  operating  points  can  be  generated  which  determines  a  receiver 
operating  characteristic  curve  or  ROC  curve.   Different  ROC 
curves  can  be  produced  by  changing  either  one  or  both  of  the 
conditional  distributions  which  implies  changing  either  the 
signal  process  or  the  noise  process. 

A  decision  rule  which  results  from  using  a  likelihood  ratio 
criterion  in  a  model  in  which  the  input  process  is  determined  by 
a  set  of  m  random  variables  can  be  expressed  in  terms  of  a  set 
n  as  follows:  n  =  {  (ylf . . . ,ym)  :  L(ylf . . . ,ym)  >  K  }  where  K 
is  specified  in  the  same  way  as  in  the  case  in  which  m  =  1. 


III.   Three  Binary  Detection  Models 

Three  detection  models  are  examined  in  this  section. 
For  the  first  two  detection  models,  the  input  stochastic  process 
for  an  observation  is  defined  by  a  time  sequence  of  continuous 
random  variables.   The  random  variables  represent  a  sample  from  a 
continuous  parameter  stochastic  process  which  is  sampled  at  times 
such  that  the  random  variables  are  independent.   For  the  third 
detection  model,  the  input  stochastic  process  is  a  counting 
process  and  it  is  defined  by  a  single  discrete  random  variable 
that  is  equal  to  the  number  of  events  that  are  counted  during  the 
observation. 

Model  I:   In  the  first  detection  model,  a  sampled  noise 
value  is  a  value  of  a  normally  distributed  random  variable  with 
mean  zero  and  with  known  variance  a2      and  a  sampled  signal  value 
is  a  known  value  of  a  deterministic  variable.   Thus,  the  input 
process  corresponding  to  an  observation  consists  of  some  number 
m  of  independent  normal  random  variables  Y^, . . . ,Ym     each  with 
variance  a2  .   And,  for  i  =  1,2,3,  •••  ,  m,  when  a  signal  is  not 
present  the  mean  of  Y^   is  zero  and  when  a  signal  is  present  the 
mean  is  s^.   The  result  of  using  a  likelihood  ratio  decision 
rule  in  the  model  can  be  expressed  in  terms  of  a  random  variable 
Z.   This  random  variable  is  called  a  crosscorrelation  statistic 
and  it  is  defined  by   Z  =  2  Sj/Y^  where  the  sum  index   1  =  1,2, 
•••  ,m  here  and  in  the  remainder  of  this  section.   However,  it 
is  more  convenient  to  express  the  result  in  terms  of  a  statistic 
V  which  is  defined  by  V  =  Z/az.      In  terms  of  this  random 


variable,  the  conditional  probabilities  pf  and  p<j  are  given 
by:   pf  =  1  -  #(v*)   and  p<j  =  1  -  #(v*  -  dH)      where  # 
symbolizes  the  standard  normal  cumulative  distribution  function, 
the  threshold  v*  =  (l/az) • (a2 In  K  +  (1/2)  2  s^)  and 
d  =  S  Sj/a2  .   The  parameter  d  is  called  the  detection  index. 

Often,  the  input  stochastic  process  represents  a  quantity 
whose  square  is  proportional  to  power.   In  such  a  case,  the 
average  receiver  input  power  is  the  random  variable  S  Yj/m.   If 
a  signal  is  not  present,  the  expected  average  receiver  input 
power  is  N  =  S  o1  /m  =  a2      where  N  is  called  the  noise.   The 
average  receiver  input  signal  power  is  S  =  E  Sj/m  where  S  is 
called  the  signal.   In  these  terms,   d  =  m* (S/N)  where  S/N  is 
called  the  signal-to-noise  ratio. 

If  a  receiver's  input  can  be  considered  to  be  a  time 
sequence  of  continuous  voltage  values  such  as  in  the  case  of  a 
sonar  receiver,  in  some  cases  a  frequency  representation  can  be 
used  that  involves  the  concept  of  receiver  bandwidth.    In  these 
cases,  the  noise  process  is  assumed  to  be  such  that  m  =  t/<St 
where  t   is  the  integration  time   (the  duration  of  an 
observation)   and  St     is  the  time  between  samples  with 
St   =  1/[2(BW)]   where  BW  is  the  bandwidth  and  St      is 
determined  by  the  sampling  theorem.   This  implies  that  the 
detection  index  can  be  written  as  d  =  2t» (BW) • (S/N) .   By 
defining  N0  as  the  power  spectral  density  where  Nq  =  N/BW   the 
detection  index  can  also  be  written  as  d  =  2t- (S/N0) •   In 
Reference  2,  the  conditions  required  by  the  first  model  are 
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called  Case  I   and  in  the  following  sections  the  first  model  is 
called  the  Case  I  model.   A  receiver  that  processes  data  such 
that  it  would  implement  a  likelihood  ratio  decision  rule  under 
the  conditions  of  the  first  model  is  called  a  crosscorrelation 
detector.    If  the  description  of  the  input  noise  is  adequate,  a 
Case  I  model  can  be  used  to  obtain  an  estimate  of  an  upper  bound 
on  a  detection  system's  performance,  since  all  the  information 
necessary  to  define  the  signal  is  assumed  to  be  known. 

Model  II:   In  the  second  detection  model,  a  sampled  noise 
value  is  a  value  of  an  independent  normal  random  variable  with 
mean  zero  and  with  variance  a2  .   And  a  sampled  signal  value  is 
an  independent  random  variable  with  mean  zero  but  with  variance 
as.   Thus,  the  input  process  corresponding  to  an  observation 
consists  of  some  number  m  of  independent  normal  random 
variables  Yi/««»»Ym  each  with  mean  zero  and  each  with  variance 
a2      when  a  signal  is  not  present  and  with  variance  a2    +   as  when 
a  signal  is  present.   The  result  of  applying  a  likelihood  ratio 
decision  rule  in  this  model  can  be  expressed  in  terms  of  a 
statistic  X  which  is  defined  by  X  =  Z  Y^.   When  a  signal  is 
not  present,  the  statistic  X/N  has  a  chi-square  distribution 
with  m  degrees  of  freedom.   When  a  signal  is  present,  the 
statistic  X/(N+S)   has  a  chi-square  distribution  with  m 
degrees  of  freedom.   So,  in  terms  of  these  statistics,  the 
conditional  probabilities  pf  and  p^  are  defined  as  follows: 
Pf=  p(xm  *  X*/N)   and  pd  =  P{X^  >  (x*/N) • [1/ (1+S/N) ] } 


where  X^  is  a  chi-square  random  variable  with  m     degrees  of 
freedom,  x*  is  a  number  which  is  determined  by  the  decision  rule 
and  S/N  is  the  signal-to-noise  ratio.   A  receiver  that  would 
implement  a  likelihood  ratio  decision  rule  under  the  conditions 
of  the  second  model  is  called  an  energy  detector  (or  square  law 
detector) . 

The  mean  of  a  chi-square  random  variable  with  m  degrees  of 
freedom  is  m  and  the  variance  is  2m.   By  the  central  limit 
theorem,  as  the  number  of  degrees  of  freedom  of  a  chi-square 
random  variable  becomes  large  it  can  be  approximated  by  a  normal 
random  variable  with  the  same  mean  and  variance.   Hence,  for  a 
sufficiently  large  sample  size  m  or  equivalently,  for  a 
sufficiently  large  time  bandwidth  product  t* (BW) ,   the 
conditional  probabilities  pf  and  p^  can  be  approximated  by 
pf  =  1  -  *(v*)   and  pd  =  1  -  #{ [1/(1+S/N) ] • (v*  -  d2*)  }   where  the 
threshold  value  v*  =  (x*  -  ma2 )/[ (2m) ^a2 ]   and  d  =  t»(BW)(S/N)2 
if  the  concept  of  bandwidth  is  applicable.   If,  in  addition,  the 
noise  N  is  significantly  larger  than  the  signal  S,  that  is,  if 
t« (BW)  »  1  and  S/N  «  1,   then  pf  and  pd  can  be 
approximated  by:   pf  =  1  -  *(v*)   and  p<|  =  1  -  *(v*  -  d"*) . 
In  Reference  2,  the  conditions  required  for  this  approximation 
are  called  Case  II  and  in  the  following  sections  the  limiting 
form  of  the  second  model  is  called  the  Case  II  model.   In 
Reference  2,  an  example  is  given  that  uses  the  Case  II  model  to 
estimate  the  performance  of  a  passive  sonar  system. 
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Model  III:   In  the  third  detection  model,  a  sampled  noise 
value  and  a  sampled  signal  value  are  values  of  independent  random 
variables  that  are  determined  by  independent  Poisson  processes 
observed  for  a  time  interval  t.   The  noise  process  is 
characterized  by  a  counting  rate  a,  the  signal  process  is 
characterized  by  a  counting  rate  as  and  the  noise  and  signal 
processes  are  additive.   Consequently,  when  the  input  is  noise 
alone,  the  input  is  the  value  of  a  random  variable  Y  that  has  a 
Poisson  distribution  with  parameter  a»t.   And,  when  a  target  is 
not  present,  the  input  is  a  value  of  a  random  variable  that  has  a 
Poisson  distribution  with  parameter   (a  +  as) »t.   For  a 
likelihood  ratio  decision  rule,   pf  =  1  -  P(y*;a«t)   and 
Pd  =  1  -  P[y*;(ot  +  as)  *t]   where  P(y;8)   symbolizes  the 
cumulative  Poisson  distribution  function  with  parameter  9  and 
y*  is  a  threshold  value  of  y  determined  by  the  decision  rule. 
In  terms  of  the  above  notation,  both  the  mean  and  variance  of  a 
Poisson  random  variable  are  equal  to  8.   When  8  is  large,  the 
cumulative  Poisson  distribution  function  can  be  approximated  by 
the  cumulative  distribution  function  of  a  normal  random  variable 
having  the  same  mean  and  variance.   Using  this  approximation  for 
cases  where  a*t  is  large  gives   pf  =  1  -  *(v*)   and 
pd  =  1  -  *{[1/(1  +  Os/a)^](v*  -  d^)}   where  v*  =  (y*  -  at) /(at)3* 
and  d  =  a»t»(as/a)2.   In  addition,  if  the  signal  counting  rate 
is  significantly  smaller  than  the  noise  counting  rate,  that  is, 
if  a»t  »  1  and  ag/o  «  1,   then  pf  and  p<j  can  be 
approximated  by:   pf  =  1  -  •(v*)   and  Pa  =  1  -  *(v*  -  d^) . 
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This  detection  model  might  be  used  to  describe  a  receiver  whose 
input  for  an  observation  is  the  number  of  photons  counted  by  a 
radiation  detector  in  situations  where  act,      the  expected  number 
of  counts  when  no  signal  is  present,  is  of  the  order  of  thirty  or 
more.   In  this  case,  the  target  might  be  a  wake  segment  for 
example. 

When  a  likelihood  ratio  decision  rule  is  used  in  the  three 
models  discussed  above,  for  the  first  model  and  under  limiting 
conditions  for  the  second  and  third  models,  the  following  result 
is  obtained:   pf  =  1  -  #(v*)   and  p<j  =  1  -  #(v*  -  d1*)      where  the 
definition  of  v*  depends  on  the  noise  power  N  for  the  first 
and  second  models.   For  a  sonar  receiver  described  by  the  first 
model,  that  is,  by  the  Case  I  model:   d  =  2t» (BW) (S/N) .   For  a 
sonar  receiver  described  under  the  limiting  conditions  for  the 
second  model,  that  is,  by  the  Case  II  model,   d  =  t» (BW) (S/N)2 . 
So,  in  either  a  Case  I  model  or  a  Case  II  model  of  a  sonar 
receiver,  the  detection  index  d  is  a  function  of  the  time 
bandwidth  product  t» (BW)   and  signal-to-noise  ratio  S/N.   Since 
sonar  equations  relate  S/N  to  system,  target  and  environmental 
parameters,  a  sonar  equation  can  be  used  to  relate  S/N  to  these 
parameters  in  a  model  of  a  sonar  receiver. 
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IV.  General  Detection  Models 

The  detection  models  that  have  been  considered  to  this 
point  are  based  on  binary  detection  theory.   After  each 
observation,  a  receiver  decides  either  that  the  input 
corresponding  to  the  observation  was  noise  or  else  it  decides  it 
was  signal  plus  noise.   However,  in  some  detection  systems  this 
decision  is  delayed.   In  a  computational  sense,  a  model  of  such  a 
detection  system  is  generally  complex  relative  to  a  binary 
detection  model.   To  illustrate  this,  consider  an  active  sonar 
system  whose  receiver  includes  an  operator.   Suppose  the 
probability  that  the  operator  will  detect  a  target  echo  has  been 
determined  in  a  laboratory  experiment  in  which  the  operator  was 
required  to  decide  after  each  input  corresponding  to  a  resolution 
cell  that  either  the  input  was  a  target  echo  (signal)  and  noise 
or  the  input  was  noise  alone.   In  addition,  suppose  that  under 
operational  conditions  the  operator  normally  delays  this 
decision.   Then,  in  general,  the  probability  that  the  operator 
will  decide  that  the  input  corresponding  to  a  resolution  cell 
that  contains  a  target  is  a  target  echo  and  noise  will  not  be 
equal  to  the  probability  of  the  event  in  the  forced  choice 
experiment.   And,  in  addition,  the  probability  that  the  operator 
will  decide  the  input  corresponding  to  a  resolution  cell  that 
does  not  contain  a  target  is  a  target  echo  and  noise  will  not  be 
equal  to  the  probability  of  this  event  in  the  forced  choice 
experiment.   Consequently,  in  general,  the  value  of  both  p^  and 
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Pf  for  an  operational  environment  will  be  different  than  that 
for  the  laboratory  environment. 

One  model  that  has  been  proposed  to  deal  with  this  kind  of 
situation  defines  the  event  that  a  receiver  decides  that  the 
input  corresponding  to  a  resolution  cell  is  signal  and  noise  to 
be  equivalent  to  the  event  that  out  of  n  consecutive 
observations  at  least  k  of  them  would  result  in  the  decision 
that  the  input  was  signal  and  noise  in  a  forced  choice 
experiment.   The  model  is  said  to  be  based  on  an  k-out-of-n 
detection  criterion.   With  this  criterion,  the  probability  that  a 
target  will  be  first  detected  on  the  j^n  observation  can  be 
found  as  follows:   Determine  the  2J   sequences  of  forced  choice 
responses  that  could  result  for  a  sequence  of  j   consecutive 
observations.   Next,  determine  the  probability  of  occurrence  for 
each  sequence  that  first  satisfies  the  k-out-of-n  detection 
criterion  on  the  j^n  observation.   The  probability  of  first 
detection  on  the  j*-n  observation  is  equal  to  the  sum  of  these 
probabilities.   The  cumulative  probability  of  detection  at  the 
jtn  observation  is  the  sum  of  the  probabilities  of  first 
detection  on  the  ith  observation  for  i=l,2,  •••  ,  j . 
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V.  Signal -to-Noise  Ratio  Detection  Models 

In  some  radar  and  sonar  detection  models,  for  a 
specified  value  of  Pf,  a  minimum  acceptable  value  of  p^   is 
defined.   This  minimum  acceptable  value  of  p^  and  the  specified 
value  of  pf  define  what  can  be  called  a  minimum  acceptable 
signal-to-noise  ratio   (S/N)m  if  pd   is  a  nondecreasing 
function  of  S/N.   In  some  sonar  detection  models,  this 
signal-to-noise  ratio  in  decibels  is  called  the  detection 
threshold  DT.   In  symbols,   DT  =  10  log (S/N) m.   If  the  minimum 
acceptable  value  of  p<j   is  .5,  then  DT   is  usually  called  the 
recognition  differential  RD.   The  difference  between  the 
signal-to-noise  ratio  in  decibels  and  RD   (or  DT)   is  called  the 
signal  excess  SE.   In  symbols,   SE  =  10  log (S/N)  -  RD. 
One  interpretation  of  signal  excess  is  that  for  a 
localization  region  containing  a  target  detection  occurs  with 
probability  one  if  SE  >  0  and  with  probability  zero  if  SE  <  0. 
This  interpretation  provides  the  basis  for  defining  detection  in 
the  three  encounter  detection  models  that  are  discussed  in 
Section  VII.   A  more  consistent  interpretation  of  signal  excess 
is:   If  SE  >  0,   then  the  probability  of  detection  p<j  is 
greater  than  or  equal  to  the  minimum  acceptable  value   (.5   if 
recognition  differential  RD  is  used  to  define  signal  excess) . 
For  cases  where  p<j  increases  rapidly  with  signal  excess  in  the 
neighborhood  of  zero  signal  excess,  the  two  interpretations  are 
essentially  equivalent.   For  a  discussion  of  this  point  as  well 
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as  a  discussion  of  an  operational  case  in  which  receiver 
decisions  are  delayed,  see  Reference  3. 

Signal  excess  (signal-to-noise  ratio)  detection  models 
provide  a  basis  for  general  detection  models,  in  particular, 
models  that  describe  nonstationary  noise  and  signal  processes  and 
randomly  changing  decision  rules.   This  is  illustrated  by  the 
models  described  in  Section  VII.   In  addition,  signal  excess 
models  provide  a  basis  for  delayed  receiver  decision  models. 
This  is  illustrated  by  the  active  sonar  detection  models  in  both 
Reference  4  and  Reference  5  that  are  based  on  a  k-out-of-n 
detection  criterion.   In  all  of  these  models,  the  signal-to-noise 
ratio  and  the  recognition  differential  are  random  variables. 

Using  X(t)   to  represent  a  random  variable  corresponding  to 
an  index  time  t  and  a  subscript  to  identify  the  random  variable 
in  such  models,  for  a  passive  sonar  receiver,  the  signal-to-noise 
ratio  in  decibels  associated  with  a  decision  at  the  index  time 
is:   XSL(t)  -  XTL(t)  -  [XNL(t)  -  XDI(t)].   In  this  expression, 
SL  represents  source  level,   TL  represents  transmission  loss, 
NL  represents  noise  level  and  DI  represents  directivity  index. 
Since  signal  excess  SE  is  defined  to  be  the  difference  in 
decibels  between  the  signal-to-noise  ratio  and  the  recognition 
differential  (or  detection  threshold) ,  it  too  is  a  random 
variable  and,  for  any  decision  time  t,  one  can  write: 
(1)   XSE(t)  =  XSL(t)  -  XTL(t)  -  [XNL(t)  -  XDI(t)]  -  XRD(t). 
The  distributions  of  the  random  variables  on  the  right  side  of 
Equation  1  determine  the  distribution  of  the  signal  excess.   In 
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the  passive  sonar  detection  model  described  in  Reference  6, 
xSL(fc) '   xRD(t)   and,  in  effect,   XNL(t)   are  normally 
distributed  random  variables  while  X<jL(t)   is  a  uniformly 
distributed  random  variable.   In  the  three  signal  excess  models 
that  are  described  in  Section  VII,  all  of  the  random  variables  in 
Equation  1  are  normally  distributed. 

It  is  sometimes  convenient  to  write  Equation  1  as  follows: 

(2)  XSE(t)  =  SE(t)  +  X(t). 

In  Equation  2,   SE(t)   is  the  expected  value  of  the  signal  excess 
determined  by  the  following  expected  value  equation: 

(3)  SE(t)  =  SL(t)  -  TL(t)  -  [NL(t)  -  DI(t)]  -  RD(t) 

where  each  term  on  the  right  represents  the  expected  value  of  the 
indicated  random  variable  and  X(t)   is  a  random  variable  that 
determines  the  stochastic  character  of  the  signal  excess.   Since 
SE(t)   is  the  mean  of  XSE(t) ;  by  Equation  2,  the  mean  of  X(t) 
is  equal  to  zero  and  the  standard  deviation  of  X(t)   is  equal  to 
the  standard  deviation  of  XSE(t) .   If  a     represents  the 
standard  deviation  of  XsE(t)   and  tne  random  variables  on  the 
right  side  of  Equation  1  are  statistically  independent,  then 
°2    =  aSL  +  aTL  +  aNL  +  aDI  +  aRD*   This  relation  has  been  used  to 
determine  a  standard  deviation  for  the  signal  excess  in 
operational  models. 
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VI.  General  Encounter  Models 

A  basic  problem  associated  with  search  modeling  is  that 
of  determining  the  probability  that  a  target  will  be  detected  by 
a  detection  system  during  an  encounter  with  one  or  more  detection 
systems.   In  the  encounter  models  that  are  considered  in  this 
report,  during  a  search,  observations  are  made  of  a  series  of 
localization  regions.   For  a  localization  region  that  contains  a 
target,   p<j   is  the  probability  of  the  event  {D±   n  H^},  that  is 
the  event  detection.  For  a  localization  region  that  does  not 
contain  a  target,  pf   is  the  probability  of  the  event   {D^  n  Hq>, 
that  is,  the  event  false  alarm.   In  addition,  in  these  models, 
the  time  to  resolve  a  false  alarm  is  ignored.   However,   p^  and 
Pf  are  assumed  to  be  determined  by  some  criterion  such  that  pf 
is  less  than  one. 

Using  the  order  number  of  a  decision  rather  than  its  time  as 
an  index  relative  to  detection  decisions  for  localization  regions 
that  contains  a  target  and  a  random  variable  N  to  represent  the 
decision  order  number  at  which  detection  first  occurs,   the 
probability  of  detection  during  an  encounter  can  be  written  as: 
P(N  <  n)  =  P(N  <  m)  +  P(N  =  m+1)  +  •••  +  P(N  =  n)   or  equally  as 
P(N  <  n)  =  1  -  [1  -  P(N  <  m)]-(l  -  qm+1)    •••  (1  -  gn)   where 


g^  =  P(N  =  i|N  <  i-1)   is  the  probability  of  the  event  detection 
at  the  itn  decision  conditioned  on  the  event  no. detection  at  an 
earlier  decision  and   1  <  m  <  n. 

The  second  expression  above  for  P(N  <  n)   is  generally  of 
greater  interest  than  the  first  expression,  since  g^  can 
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usually  be  more  directly  related  to  operational  parameters  such 
as  target  range  and  target  environment  that  determine  a  target's 
detectability  than  can  P(N  =  i) .  If  a  target  is  not  present  in 
the  localization  region  that  corresponds  to  a  decision,  then  the 
event  {D^  n  H^}  can  not  occur  and  the  probability  of  this  event 
will  be  zero. 

With  a  time  rather  than  the  order  number  to  index  a  decision 
and  a  random  variable  T  to  represent  the  time  index  at  which 
detection  first  occurs,   P(N  <  n)   becomes  P(T  <  tn)   with 
P(T  <  tn)  =  1  -  [1  -  P(T  <  tm)].[l  -  g(tm+1)]  -.  [1  -  g(tn)] 


where  now  g^  =  P(T  =  tjjT  Stj-i)  . 

If  g(tjj  «  1  for  1  =  1,2,  ..«  ,n,   then,  to  a  first 
approximation,   ln[l  -  <j(tj_)  ]  =  -g(t^)   for  1  =  1,2,  •••  ,n  and 
P(T  <  tn)  =  1  -  [1  -  P(T  <  tm)]«exp[-S  g(ti)].   This  follows 
since   P(T  <  tn)  =  1  -  [1  -  P(T  <  tm)]«exp[E  ln[l  -  g(t^) ]   where 
the  sum  index  i  =  m+1,  •••  ,n.   A  continuous  analog  to  this 
approximation  can  be  used  to  describe  an  encounter  for  which 
g(tjj  «  1  for  i  =  m,m+l,  •••  ,n  and  decisions  during  the 
encounter  can  be  considered  to  occur  continuously.   That  is,  the 
time  of  an  observation  corresponding  to  a  decision  and  the  time 
between  decisions  are  both  negligible  relative  to  the  time  of  the 
encounter. 

The  analog  can  be  developed  as  follows:   First,  let  6t  be 
the  time  between  decisions,  then  t^  =  i«6t  and  the  probability 
of  detection  P(T  <  tn)  =  1  -  [1  -  P(T  <  tm) ] «exp[-S  t (t^  » St] 
where  T(t^)  =  (l/6t) *g(ti)   is  a  detection  rate  function 
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(a  probability  of  detection  per  unit  time)   and,  in  terms  of  6t, 


the  probability   g(ti)  =  P[T  =  i»*t|T  <  (i-l)«St]. 

If  T  is  considered  to  be  a  continuous  random  variable,  the 
expression  for  P(T  <  tjj   above  indicates  that  the  sum  in  the 
exponent  should  be  replaced  by  an  integral  whose  integrand  is  a 
continuous  function  r(t).   If  t (t)   can  be  determined,  then, 
with  g(t^)   as  a  guide,  the  cumulative  probability  of  detection 
P(T  <  t)   can  be  defined  by: 


(4)  r(t)  =  lim  {(l/5t)»P(t  <  T  <  t+6t|T  <  t) } 
where  the  limit  is  for  St     approaching  zero.   Equation  4 
implies  the  differential  equation:   dp(t)/dt  =  [1  -  p(t)]»r(t) 
where  p(t)  =  P(T  <  t) .   The  solution  to  this  differential 

equation  is: 

ptn 

(5)  P(T  <  tn)  =  1  -  P(T  <  tm).exp[-Jtm  r(t)dt] 

where  t  is  the  time  index  for  a  decision  during  an  encounter,   tm 
is  some  time  during  the  encounter  and  tn  >  tm.   A  r (t)   that  is 
based  on  a  visual  detection  model  is  described  in  Reference  7. 
If  the  detection  capability  of  a  detection  system  is  assumed  to 
depend  on  a  target's  position  relative  to  the  detection  system 
during  an  encounter  but  not  to  depend  on  the  clock  time,  then  the 
time  index  of  a  decision  can  be  a  relative  index  that  determines 
the  target  position  that  is  associated  with  a  decision  rather 
than  the  clock  time  associated  with  the  decision. 

The  above  results  apply  to  the  case  of  an  encounter  between 
a  target  and  a  collection  of  detection  systems.  However,  if  the 
detection  systems  are  not  collocated,  it  is  generally  convenient 
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to  describe  encounters  of  this  kind  in  terms  of  encounters 
between  the  target  and  the  individual  detection  systems.   In 
either  case,  if  the  event  target  detection  for  a  detection  system 
is  not  independent  of  the  event  for  other  detection  systems,  then 
in  order  to  describe  this  in  an  encounter  model  the  correlation 
between  the  input  to  the  detection  system  and  the  inputs  to  the 
other  detection  systems  must  be  specified.   This  has  been  done  in 
some  models  as  follows:   First  determine  the  probability  of 
detection  for  each  system  acting  alone.   Let  P^  be  the 
probability  that  the  ith  system  detects  the  target  during  the 
encounter  under  this  condition.   Next,  consider  two  cases:   In 
the  first  case,  the  random  factors  that  determine  detection  for  a 
system  are  independent  of  those  that  determine  detection  for  the 
remaining  systems.   In  the  second  case,  the  random  factors  that 
determine  detection  for  the  systems  are  completely  dependent.   In 
the  first  case,  the  probability  that  at  least  one  system  detects 
the  target  is  given  by:   Pj  =  1  -  (1  -  Px)  •  (1  -  P2)  •••  (1  -  Pn) 
where  n   is  the  number  of  detection  systems  involved.   In  the 
second  case,  the  probability  that  none  of  the  systems  detect  the 
target  is  given  by:  1  -  PD  -  1  -  Pm  where  Pm  >  P^   for 
i  =  1,2,  •••  ,n  since  if  the  m^h  system  does  not  detect  the 
target,  none  of  the  remaining  systems  will  detect  it.   The 
probability  that  at  least  one  system  detects  the  target  is  given 
by:   P  =  OfP[)  +  (1  -  crJ'Pj  where  a  determines  the  degree  of 
correlation  and  0  <  a  <   1.   A  way  to  determine  a  value  for  a 
is  described  in  Reference  8. 
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VII.   Three  Signal  Excess  Encounter  Models 

In  the  three  models  described  in  this  section, 
detection  is  defined  in  terms  of  signal  excess  as  it  is  in 
Section  V.   Each  of  the  models  determines  a  cumulative 
probability  of  detection  for  a  target  in  an  encounter  with  a 
passive  sonar  system.   An  observation  in  the  models  is  indexed  by 
time  and  the  index  can  usually  be  considered  to  be  the  time  at 
the  end  of  the  observation.   During  an  encounter,  observations 
are  made  of  a  series  of  localization  regions.   For  localization 
regions  that  do  not  contain  a  target,  the  signal  observed  during 
the  observation  of  the  region  is  zero.   The  time  to  resolve  a 
false  alarm  is  zero.   However,  false  alarms  are  not  ignored  in 
that  the  value  of  RD   (or  DT)   is  determined  by  some  specified 
false  alarm  probability.   Consequently,  although  the  time  to 
resolve  a  false  alarm  is  zero,  the  cost  associated  with  a  false 
alarm  is  not  zero.   In  determining  signal  excess  in  the  models, 
it  is  convenient  to  use  Equation  2.   For  each  decision  in  an 
encounter,  there  is  a  random  variable  X(t)   defined  by  Equation 
2  that  determines  the  random  character  of  the  signal  excess. 
For  a  sequence  of  decisions,  the  set  of  these  random  variables 
ordered  by  their  time  index  constitutes  a  stochastic  process. 
The  joint  distributions  of  these  random  variables  determines  the 
nature  of  this  process.   In  the  three  encounter  models  described 
in  this  section,  the  stochastic  process  is  called  a  lambda-sigma 
jump  process.   The  jumps  in  the  series  occur  at  times  determined 
by  a  Poisson  process  with  a  mean  rate  lambda.   This  implies  that 
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the  time  between  jumps  is  a  random  variable  with  an  exponential 
distribution  and  that  the  expected  times  between  jumps   r   is 
equal  to  the  reciprocal  of  lambda.   The  time  series  that  are 
generated  by  lambda-sigma  jump  processes  are  represented  by 
Figure  2   below. 


dB 


time 


Figure  2.  A  time  series  representing  a  realization  of  a  lambda- 
sigma  jump  process.  On  the  plot,  a  in  dB  equals  one  unit  on 
the  vertical  axis  and  r  equals  one  time  unit  on  the  horizontal 
axis. 

From  Figure  2,  note  that  the  observed  values  of  neighboring 
random  variables  are  equal  unless  a  jump  has  occurred  between 
them.   When  a  jump  occurs,  the  first  random  variable  after  the 
jump  is  normally  distributed  with  mean  zero  and  variance  a2      and 
it  is  independent  of  all  the  random  variables  before  the  jump. 
Conditioned  on  a  jump  pattern,  this  random  variable  and  all  the 
random  variables  between  it  and  the  next  jump  are  dependent  and 
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the  correlation  coefficient  between  any  pair  is  one.   That  is,  if 
the  value  of  the  signal  excess  is  known  at  some  time,  then  all  of 
the  values  between  the  last  jump  before  that  time  and  the  first 
jump  after  that  time  are  also  known.  However,  since  the  jumps 
occur  randomly,  knowing  the  value  of  the  signal  excess  with 
certainty  at  some  time  does  not  determine  the  values  of  the 
signal  excess  with  certainty  at  neighboring  times.   In  the 
unconditioned  case,  the  correlation  coefficient  between  the 
random  variables  X(t)   and  X(t+T)   is  equal  to   1/e.   For  this 
reason,  r      is  referred  to  as  a  relaxation  time.   If  only  the 
transmission  loss  is  considered  to  be  determined  by  a  stochastic 
process,  then  only  the  transmission  loss  will  be  a  random 
quantity  and  t   is  the  expected  time  between  jumps  in  the  random 
component  of  the  transmission  loss.   It  appears  that  the  use  of 
the  lambda-sigma  jump  process  is  based  more  on  past  practice  than 
on  experimental  justification.   In  this  regard,  see  Reference  9. 
By  referring  to  Equation  1,  it  can  be  seen  that  the  lambda-sigma 
jump  process  is  equal  to  the  sum  of  the  stochastic  processes  that 
determine  the  random  variables  on  the  right  side  of  this 
equation.   But,  the  nature  of  the  processes  should  depend  on  the 
encounter  that  is  being  modeled.   This  suggests  that  an  adequate 
description  of  an  encounter  using  a  lambda-sigma  jump  process 
might  not  be  possible  for  some  encounters. 

In  the  three  encounter  models  described  below,  detection  is 
defined  in  terms  of  signal  excess  as  described  in  Section  IV  and 
decisions  are  indexed  by  a  time  that  can  usually  be  considered  to 
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be  the  time  of  the  decision.   During  an  encounter,  observations 
are  made  of  a  series  of  localization  regions.   For  localization 
regions  that  do  not  contain  a  target,  the  signal  observed  during 
an  observation  of  the  region  is  zero.   In  the  models,  the  time  to 
resolve  a  false  alarm  is  zero.   However,  false  alarms  are  not 
ignored  in  that  the  value  of  RD   (or  DT)   is  determined  by  some 
specified  false  alarm  probability.   Consequently,  although  the 
time  to  resolve  a  false  alarm  is  zero,  the  cost  associated  with  a 
false  alarm  is  not  zero. 

The  First  Passive  Sonar  Encounter  Detection  Model:   This 
model  describes  an  encounter  in  terms  of  a  series  of  decisions 
with  each  decision  based  on  the  signal  excess  XSE(t)   at  a  time 
corresponding  to  the  end  of  an  observation.   The  observations  are 
of  equal  duration  and  the  integration  time  that  determines  the 
recognition  differential  is  equal  to  the  duration  of  the 
observations.   In  the  model,   XsE(t)   i-s  determined  by  a  lambda- 
sigma  jump  process  .    For  an  encounter  in  which  SE(t)   is 
unimodal  and  in  which  the  time  of  the  single  maximum  is  prior  to 
the  end  of  the  encounter,  it  is  shown  in  Reference  10  that  the 
probability  p  that  detection  will  occur  during  the  encounter 
is  given  by  the  following  equation: 

(6)    p  =  1  -  [(1  -  pc)/(l  -  /9-pc)]-(l  -  /3-Px)  •••  (1  -  /3-Pm) 
where  /3  =  1  -  exp(-<St/r)   and  p^  =  *[SE(ti)/a]   for  i  =  1,2, 
•••  ,m.   Here,  St      indicates  the  duration  of  an  observation  and 
*   indicates  the  standard  normal  cumulative  distribution  function 
as  before.   The  integer  c  is  the  index  of  the  decision  time  tc 
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for  which  SE(tc)   is  greater  than  or  equal  to  SE(ti)   for  any 
time  t£  and  t^  <  tc  £  tm. 

As  r  approaches  zero,   /9  approaches  one  and  Equation  6 
approaches  this  form: 

(7)  p  =  1  -  (1  -  Pi)  ...  (1  -  i^). 

In  this  limit,  the  signal  excess  random  variables  are  all 
independent.   Note  that  Equation  7  applies  without  the  condition 
that  SE(t)   be  unimodal. 

As  t  approaches  infinity,   /3  approaches  zero  and 
Equation  6  approaches  this  form: 

(8)  p  =  pc. 

In  this  limit,  the  correlation  coefficient  between  any  pair  of 
signal  excess  random  variables  is  equal  to  one.   Note  that 
Equation  8   applies  without  the  condition  that  SE(t)   be 
unimodal.   Equation  8  defines  a  complete  dependence  encounter 
model . 

The  Second  Passive  Sonar  Encounter  Detection  Model:   This 
model  is  in  a  sense  a  third  limiting  form  of  the  first  passive 
sonar  encounter  detection  model.   In  this  limit,  the  time  between 
decisions  approaches  zero.   However,  in  this  limit  the 
integration  time  that  determines  the  recognition  differential  is 
not  equal  to  St     and  it  does  not  approach  zero.   It  is,  in 
effect,  chosen  by  the  user  of  the  model  through  the  user's  choice 
of  the  value  for  the  recognition  differential.    For  an  encounter 
that  begins  at  tl  and  ends  at  t2  and  for  which  XSE(t)   is 
determined  by  a  lambda-sigma  jump  process  and  SE(t)   is 
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unimodal,  it  is  shown  in  Reference  10  that  for  this  limit, 
Equation  6  has  the  following  form: 

rt2 

(9)  p  =  1  -  [1  -  p(tc)].exp[-(l/T).Jtl  p(t)  dt] 

where  p(t)  =  *[SE(t)/a]   and  where  now  tc   is  the  encounter 
time  such  that  SE(tc)   is  greater  than  equal  to  SE(t)   for  any 
other  encounter  time  t  and   tl  <  tc  <  t2. 

The  Third  Passive  Sonar  Encounter  Detection  Model:    This 
model  describes  an  encounter  between  a  target  and  a  passive  sonar 
detection  system  in  which  detection  occurs  during  an  encounter  if 
the  average  value  of  the  square  of  the  continuously  observed 
signal-to-noise  ratio  over  a  time  interval  of  length  u   is 
greater  than  or  equal  to  the  square  of  the  signal-to-noise  ratio 
that  determines  the  recognition  differential  for  an  integration 
time  equal  to  u.   With  R(s)   the  random  signal-to-noise  ratio 
at  a  time  s  and  Rm(u)   the  random  signal-to-noise  ratio  that 
determines  the  random  recognition  differential  for  an  integration 
time  u,  detection  during  an  encounter  occurs  at  the  first  time 
t  that  the  following  inequality  is  satisfied: 

(10)  (1/u)  Jt-u  [R(s)/Rm(u)r  ds  >  1 

where  the  time  origin  is  chosen  so  that  t  >  0   and  where  the 
integration  time  u  =  t  for  t  <  t0  and  u  =  t0   for  t  >  t0 
where  t0   is  a  maximum  integration  time.   The  random  integrand 
in  the  inequality  is  related  to  the  random  signal  excess  at  the 
time  s   for  an  integration  time  u.   The  relation  is: 

(11)  10  log  [R(s)/Rm(u)]2  =  2[SE(s;u)  +  X(s)] 
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where  SE(s;u)   is  the  expected  value  of  the  signal  excess  at  a 
time  s  for  an  integration  time  u  and  X(s)   is  the  random 
component  of  the  signal  excess  at  the  time  s.   In  the  model, 
X(s)   is  determined  by  a  lambda-sigma  jump  process  and  SE(s;u) 
is  determined  by  an  expected  value  sonar  equation  with  a 
recognition  differential  RD(u)  =  10  log  rm(u)  .   Here,   rmft)   ^s 
the  value  of  the  signal-to-noise  ratio  that  gives  a  probability 
of  detection  equal  to   .5  for  an  integration  time  t  and  a 
specified  probability  of  false  alarm  pf.   With  the  signal 
detection  process  described  by  a  Case  II  signal  detection  model, 
the  detection  index  necessary  to  give  the  required  operating 
point   (pf,.5)   is  related  to  the  integration  time  t  and  the 
signal-to-noise  ratio  rm(t)   by: 

(12)  d  =  u.(BW)[rm(t)]' 

where  BW  is  the  bandwidth  of  the  receiver.   For  a  spectrum 
analyzer,   BW  would  be  the  bandwidth  corresponding  to  a  given 
frequency  resolution  and  d  would  be  the  detection  index 
required  in  order  to  be  at  the  operating  point   (pf,.5)   for  a 
signal  that  was  contained  within  a  bandwidth  BW.   Since  d  in 
Equation  12  must  be  the  same  for  t  =  u  and  t  =  t0, 

(13)  RD(u)  =  5  log(t0/u)  +  RD(t0) 

where  t0  is  the  maximum  integration  time.   Then,  since 
SE(s;u)  -  SE(s;t0)  =  RD(t0)  -  RD(u) ,  by  using  Equation  13  and 
Equation  11,  Relation  10  becomes: 

I*      (1/5) [X(s)  +  SE(s;t0)  -  5  log(t0)] 

(14)  Jt_u  10  ds  >  1 
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where  as  above  the  time  origin  is  chosen  so  that  t  >  0,  the 
integration  time  u  =  t  for  t  <  t0  and  u  =  t0  for  t  >  t0 
and  where  SE(s;t0)   is  the  expected  value  of  the  signal  excess 
at  the  time  s   for  a  recognition  differential  determined  by  an 
integration  time  t0.   In  an  encounter,  detection  occurs  the 
first  time  that  Relation  14   is  satisfied. 

As  is  pointed  out  in  Reference  11,  the  appeal  of  the  Third 
Passive  Sonar  Encounter  Detection  Model  relative  to  the  Second 
and  First  Passive  Sonar  Encounter  Detection  Models  is  that  it 
appears  to  more  closely  describe  the  detection  process  in  passive 
sonar  detection  systems  that  display  their  processed  data  to  an 
operator  in  a  continuous  manner  over  a  time  window  of  duration 
t0.   However,  results  reported  in  Reference  12   indicate  that  the 
difference  between  the  three  models  may  not  be  significant  in 
some  types  of  encounters. 
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VIII.   Straight  Line  Encounters 

Suppose  a  target's  detectability  depends  on  its  range 
from  a  detection  system  and  that  the  probability  of  detection  is 
effectively  zero  beyond  a  range  rm  for  any  target  azimuth.   In 
this  report,  an  encounter  between  the  target  and  the  detection 
system  is  the  event  that  the  range  between  the  target  and  the 
detection  system  is  less  than  or  equal  to  rm.   In  addition, 
suppose  rm  is  small  enough  so  that  when  the  target  and  the 
detection  system  are  having  an  encounter  they  can  be  considered 
to  be  moving  on  planes  parallel  to  a  tangent  plane  to  the  earth's 
surface  at  some  point  in  their  vicinity.   If  this  is  the  case, 
then  while  the  target  and  detection  system  maintain  a  constant 
course  and  speed  during  an  encounter,  the  encounter  is  called  a 
straight  line  encounter. 

A  straight  line  encounter  can  be  described  in  terms  of  a  two 
dimensional  rectangular  coordinate  system  whose  plane  is  parallel 
to  the  tangent  plane  to  the  earth.   If  the  coordinate  system  is 
stationary  relative  to  the  detection  system  with  the  detection 
system  located  at  the  origin  and  is  oriented  so  that  the  target's 
motion  is  parallel  to  the  y-axis  and  is  in  the  positive 
y-direction,  then  the  target's  x-coordinate  during  a  straight 
line  encounter  will  be  constant.   The  constant  is  equal  to  the 
target's  horizontal  range  at  the  closest  point  of  approach  (CPA) 
on  the  straight  line  track  on  which  the  target  is  moving  relative 
to  the  detection  system  during  the  encounter.   This  range  is 
called  the  target's  lateral  range. 
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A  complete  straight  line  encounter  is  a  straight  line 
encounter  that  begins  at  a  range  from  a  detection  system  that  is 
greater  than  or  equal  to  rm  and  continues  past  CPA  to  a  range 
from  the  detection  system  that  is  again  equal  to  or  greater  than 
rm.   Let  p(x)   be  the  cumulative  probability  that  a  target  is 
detected  by  a  detection  system  in  a  complete  straight  line 
encounter  in  which  the  target's  lateral  range  is  x.   Then  the 
function  p(x)   defines  what  is  called  a  lateral  range  curve  or 
lateral  range  function. 

Let  p  be  the  probability  that  a  target  is  detected  during  a 
complete  straight  line  encounter.   If  the  lateral  range  of  a 
target  in  a  straight  line  encounter  is  assumed  to  be  a  continuous 
random  variable  X  with  a  uniform  distribution  with  fx(x)  =  ^Z* 
for   |x|  <  a/2  and  p(x)  =  0  for   |x|  >  a/2,  then  the 
probability  that  a  target  will  be  detected  during  a  complete 
straight  line  encounter  is  given  by: 

poo 

(15)  p  =  (1/a)  J  _«,  p(x)  dx 

where  the  limits  of  integration  can  be  used  since  the  value  of 
p(x)  is  zero  for   |x|  >  a/2.   Equation  15   suggests  a  measure  of 
a  detection  system's  capability  to  detect  a  target  in  a  straight 
line  encounter.   The  measure  W  is  called  sweep  width  and 

(16)  W  -  ,)_«>  p(x)  dx. 

With  this  definition,  Equation  16  becomes:   p  =  (l/2a) «W. 
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IX.  Two  Intermittent  Signal  Encounter  Models 

In  the  intermittent  signal  encounter  models  that  are 
described  in  this  section,  an  encounter  is  a  complete  straight 
line  encounter,  and  during  an  encounter  a  target  either  emits  a 
signal  (an  acoustic  transient)  or  its  presence  (a  visible 
submarine  mast)  is  the  cause  of  a  signal  at  various  times.   Two 
cases  are  considered:   In  the  first  case,  the  signals  occur 
periodically,  the  signals  are  of  length  St     and  the  time  between 
the  occurrence  of  signals  is  t  where  r   >   St.      In  the  second 
case,  fit  =  0  (the  signals  are  instantaneous)  and  the  signals 
occur  at  times  determined  by  a  Poisson  process  for  which  the 
expected  time  between  signals  is  equal  to  r.   In  the  model,  the 
detectability  of  a  target  signal  depends  on  a  target's  horizontal 
range  from  a  detection  system,  but  on  no  other  factors.   If  a 
signal  occurs  while  the  target  is  within  a  range  r,  it  will  be 
detected.   For  a  continuous  signal,  the  lateral  range  function  of 
a  detection  system  for  a  target  is:   p(x)  =1   for   |x|  <  r  and 
p(x)  =  0   for   |x|  >  r  where  the  horizontal  range  r  is 
determined  by  the  characteristics  of  the  detection  system  and  the 
target.   The  geometry  for  an  encounter  is  shown  in  Figure  1 
below. 

For  intermittent  signals,  the  length  of  a  target's  track 
relative  to  a  detection  system  on  which  a  signal  will  be  detected 
is  2*(r2  -  x2)1  +  wit  where  w  is  the  speed  of  the  target 
relative  to  the  detection  system.   So,  a  target's  exposure  time 
during  an  encounter  is   (2/w) • (r2  -  x2  )  **  +  St. 
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For  periodic  signals,  there  are  two  cases.   In  the  first 
case,   r  >  w(t  -  $t)/2.   In  this  case,  the  signals  result  in  the 
following  lateral  range  function: 

p(x)  =  0   for   |x|  >  r 
(17)       p(x)  =  1   for   |x|  <  {r2  -  [w(t  -  5t)/2]2}^ 

p(x)  =  [2/(wr)  ]  •  (r2  -  x2  )  %   +  St/r     otherwise 
In  the  second  case,   r  <  w (t  -  6t)/2      and  the  middle  equality 
in  Equation  17   does  not  apply. 

For  signals  that  are  instantaneous  and  whose  occurrence  is 
determined  by  a  Poisson  process,  the  signals  result  in  the 
following  lateral  range  function: 

p(x)  =  1  -  exp(-[2/(wT)]«(r2  -  x2)3*}   for   |x|  <  r 

p(x)  =  0   for   |x|  >  r. 


(18) 


Figure  3.   The  encounter  geometry  for  the  two  intermittent  signal 
models  described  here. 
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For  signals  whose  occurrence  is  determined  by  a  Poisson 
process  and  for  which  St  >   0,  signals  can  overlap.   If  this  is 
allowed,  then  Equation  18  can  be  modified  to  describe  this  case 
by  adding  St/r     to  the  term  in  the  exponent  of  Equation  18  that 
is  within  the  square  brackets.   In  particular,  note  that  this 
modified  Equation  18  can  be  approximated  by  the  bottom  equality 
in  Equation  17  when   (2/wr)(r2  -  x2  )  ^  +  St/r   «   1.   This 
implies  that  when  the  expected  time  t  between  signals  is  large 
relative  to  the  exposure  time   (2/w) (r2    -  x2  ) %  +   St,    the  periodic 
occurrence  model  and  the  random  occurrence  model  are  effectively 
equivalent. 
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X.   A  Random  Search  Model 

A  search  of  a  region  in  which  a  target  moves  on  a  track 
that  consists  of  a  number  of  straight  segments  placed  in  such  a 
way  that  in  a  limiting  sense  every  section  of  the  region  is 
equally  likely  to  be  searched  on  a  segment  is  referred  to  as  a 
random  search  in  Reference  6.   Representation  of  a  search  region 
with  the  track  segments  that  could  be  imagined  to  be  the  tracks 
of  a  random  search  are  shown  in  Figure  4 . 


Figure  4.   A  search  region  and  a  track  that  could  be  described  as 
a  random  search  track. 

Two  developments  of  a  model  to  describe  this  kind  of  search 
are  contained  in  this  section.   The  first  development  is  based  on 
the  following  conditions:   1.  A  target  is  at  a  fixed  position 
within  a  defined  search  region.   2.  A  searcher's  track  is  a 
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sequence  of  straight  line  segments  that  are  within  the  search 

region.   3.  The  searcher's  detection  system  is  such  that  while  on 

a  track  segment,  a  rectangle  is  searched  that  is  contained  within 

the  search  region,  is  of  length  equal  to  the  length  of  the  track 

segment  and  is  oriented  so  that  its  long  axis  is  parallel  to  the 

track  segment.   4.  The  probability  that  the  searcher's  detection 

system  will  detect  a  target  while  on  a  track  segment  with  a 

search  rectangle  that  does  not  contain  the  target  is  zero.   The 

probability  that  the  searcher's  detection  system  will  detect  a 

target  while  on  a  track  segment  with  a  search  rectangle  that 

contains  the  target  is  p(x)   where  x  is  the  target's  lateral 

range  for  the  track  segment  and  p(x)   is  the  lateral  range  curve 

for  a  complete  straight  line  encounter  lateral  range  x.  A 

representation  of  a  search  rectangle  is  shown  in  Figure  5  below. 

5.  The  track  segments  are  located  in  such  a  way  that  the  event 

that  the  target  is  within  the  search  rectangle  associated  with  a 

track  segment  is  independent  of  the  event  that  the  target  in  the 

search  rectangle  associated  with  any  other  track  segment.   And 

the  probability  of  the  event  is  equal  to  the  ratio  of  the  area  of 

■ 
the  search  rectangle  to  the  area  of  the  search  region  and,  given 

a  target  is  within  a  search  rectangle,  its  position  is  uniformly 

distributed  over  the  rectangle. 

Condition  4   implies  that  the  random  search  model  is  based 

on  the  concept  of  a  complete  straight  line  encounter.   The 

definition  of  an  encounter  that  is  intended  here  is  that  given  in 

Section  VI.   This  implies  that  in  the  random  search  model  the 
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time  to  resolve  a  false  alarm  is  zero.   However,  for  the  model, 
Pjj  and  pf  are  considered  to  be  determined  by  some  criterion 
such  that  pf  is  less  than  one.   Consequently,  although  the  time 
to  resolve  a  false  alarm  is  zero  in  the  model,  the  cost 
associated  with  a  false  alarm  is  not  zero.   (A  concise  model  that 
accounts  for  the  time  to  resolve  false  alarms  is  described  in 
Reference  13.)   Condition  4   also  implies  that  when  a  searcher  is 
on  a  track  segment  with  a  search  rectangle  that  contains  a 
target,  the  encounter  is  a  complete  straight  line  encounter.   And 
Condition  5,  which  can  be  considered  to  specify  a  random 
arrangement  of  the  track  segments,  implies  that  when  this  is  the 
case,  for  the  complete  straight  line  encounter,  the  target's 
lateral  range  is  a  random  variable  that  is  uniformly  distributed 
between  -b/2  and  b/2  where  b  is  width  of  the  search 
rectangle  (the  dimension  of  the  rectangle  perpendicular  to  the 
associated  track  segment) . 


Swept  Area 


Track 


Figure  5.   A  track  segment  and  its  associated  search  rectangle, 
that  could  correspond  to  a  search  with  an  aircraft  mounted 
infrared  detection  system. 
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Based  on  the  above  considerations,  the  probability  that  a 
target  will  be  detected  while  a  searcher  is  on  a  track  segment 
with  an  associated  search  rectangle  that  contains  the  target  is 
given  by: 

(19)  ]_«,  p(x)  fx(x)  dx  =  W/b 

where   fx(*)  -  V*>  for  -b/2  <  x  <  b/2   and   fx(*)  ■  0   and 
p(x)  =0  otherwise.   Note  that  the  left  side  of  Equation  19 
applies  to  any  complete  straight  line  encounter  in  which  the 
target's  lateral  range  for  the  encounter  is  considered  to  be  a 
random  variable  with  a  distribution  determined  by  the  probability 
density  function  fx(x) .   If  it  is  not  given  that  the  target  is 
within  the  search  rectangle  associated  with  a  track  segment,  then 
the  unconditional  probability  that  the  target  will  be  detected  on 
the  track  segment  is  given  by:   (W/b) • (6A/A)   where  6 A  is  the 
area  of  the  search  rectangle  associated  with  the  track  segment 
and  A  is  the  area  of  the  search  region.   With   1   the  length  of 
the  rectangle,  6k   =  b»l  and  the  probability  becomes:   (W«l)/A. 
Then,  since  the  event  that  the  target  will  be  in  the  search 
rectangle  of  a  track  segment  is  independent  of  the  event  that  it 
will  be  in  the  search  rectangle  of  any  other  track  segment,  the 
probability  p  that  a  random  search  consisting  of  m  track 
segments  will  detect  the  target  is  given  by: 

1  -  [1  -  (W'l^/AHl  -  (W»l2)/A]  •••  [1  -  (W-ln)/A]   where  li 
is  the  length  of  the   ith  track  segment.   The  probability  is 
also  given  by:   p  =  1  -  exp{Z  ln[l  -  (W'l-^/A]}   where  the  sum 
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index   i  =  1,2,  •••  ,n.   If   (W-ljj/A  «  1   for  i  =  1,2,  •••  ,n, 

then  this  expression  can  be  approximated  by: 

(20)      p  =   1  -  exp[-(W-l)/A] 

where   1  =  S  1^   is  the  track  length  of  the  search.   Equation  20 

is  known  as  the  random  search  formula. 

The  second  development  of  the  random  search  formula  is  based 
on  Equation  5  and  a  detection  rate  for  a  random  search  given  by: 
r(t)  =  W«v(t)/A.   With  this  detection  rate  and  Equation  5,  the 
random  search  formula  is  given  by: 
(20a)   P(T  <  t)  =  1  -  exp  {-[W»l (t) ]/A} 

where   l(t)   is  the  track  length  for  a  random  search  that  starts 
at  time   0   and  ends  at  time  t  and 


l(t)  =  Jo  v 


(20b)       l(t)  =  J0  v(s)ds. 
Replacing  P(T  <  t)   by  p  and  l(t)   by   1  gives  Equation  20. 
In  the  form  of  Equation  20a,  the  random  search  formula  indicates 
explicitly  the  relation  between  the  probability  of  detection  and 
the  duration  of  a  random  search.   Note  that  Equation  20a  implies 
that  the  sweep  width  is  independent  of  speed  over  the  range  of 
speeds  in  the  encounter. 

Reference  14  contains  an  example  of  an  application  of  the 
method  used  in  the  second  development  of  the  random  search 
formula  to  a  random  search  where  the  search  region  expands  with 
time. 
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XI.   Ladder  and  Barrier  Search  Models 

In  some  barrier  searches,  the  barrier  search  track  is  a 
ladder  search  track  relative  to  a  reference  system  that  moves 
with  the  target.   This  fact  is  used  in  the  barrier  search  model 
development  that  follows  the  two  ladder  search  model  developments 
below.   The  first  ladder  search  model  is  referred  to  as  an  ideal 
ladder  search  model  because  of  the  idealizations  that  are 
involved  in  its  description  of  a  ladder  search.   The  second 
ladder  search  model  is  referred  to  as  a  degraded  ladder  search. 
It  can  be  considered  to  describe  a  ladder  search  track  in  which 
navigational  errors  result  in  omissions  and  overlaps  in  coverage. 

An  Ideal  Ladder  Search  Model:   The  model  is  based  on  the 
following  conditions:   1.  A  ladder  search  region  is  a  rectangle 
that  contains  a  fixed  target.   2.  During  a  search  of  the  region, 
the  searcher's  detection  system  moves  on  a  set  of  m  parallel 
track  segments  of  length  b  separated  by  a  distance  s.   3 .  As 
the  detection  system  moves  along  a  track  segment,  it  searches  a 
rectangular  strip  of  length  b  and  width   s  within  the  search 
region.   4.  The  m  rectangular  strips  that  correspond  to  the  m 
track  segments  completely  cover  the  ladder  search  region  with  no 
overlap.   5.  If  a  target  is  within  the  rectangular  strip 
corresponding  to  a  track  segment,  then  there  will  be  a  complete 
straight  line  encounter  between  the  target  and  the  detection 
system  when  the  detection  system  moves  along  the  track  segment 
and  the  lateral  range  of  the  encounter  will  be  uniformly 
distributed  across  the  width  of  the  strip.   If  the  target  is  not 
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in  the  rectangular  strip,  then  there  will  not  be  an  encounter  and 
the  probability  that  the  target  will  be  detected  while  the 
detection  system  is  on  the  track  segment  is  zero. 


Figure  6.   A  schematic  representation  of  a  ladder  search  geometry 
for  a  case  in  which  the  ladder  search  track  segments  are 
superimposed  on  and  bisect  their  corresponding  rectangular 
strips. 

Since  targets  outside  of  the  rectangular  strip  that 
corresponds  to  a  track  segment  cannot  be  detected  while  a 
detection  system  is  on  the  track  segment  because  of  Condition  5, 
in  the  model,  the  sweep  width  W  of  a  searcher's  detection 
system  must  satisfy  the  relation  W  <  s  .   In  particular,   W  =  s 
only  holds  when  the  detection  system  detects  a  target  that  is  in 
a  rectangular  strip  with  probability  one  for  any  target  lateral 
range.   This  kind  of  detection  system  is  sometimes  referred  to  as 
a  cookie  cutter  detection  system.    However,  this  terminology  can 
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be  misleading  since  it  suggests  the  detection  system  detects 
equally  well  for  all  azimuths.   But  this  is  not  a  requirement  on 
the  system  in  order  that  W  =  s. 

The  ideal  ladder  search  model  implies  that  if  the  conditions 
of  the  model  are  satisfied,  then  the  probability  p  that  a 
target  will  be  detected  by  a  an  ideal  ladder  search  is  given  by: 

(21)  p  =  W/s 

where  W/s  <  1.   The  quantity  W/s   is  called  the  coverage 
factor  in  this  case. 

A  Degraded  Ladder  Search  Model:   The  above  model  implies 
perfect  navigation  in  addition  to  other  idealizations.   A  model 
of  a  ladder  search  is  given  in  Reference  6  that  could  be  used 
for  cases  in  which  this  is  a  poor  assumption.   The  model  which  is 
referred  to  here  as  a  degraded  ladder  search  model  can  be 
considered  to  describe  navigational  inaccuracies  in  terms  of 
omissions  and  overlaps  of  the  rectangular  strips.   It  can  be 
developed  as  follows:   Consider  a  random  search  in  the  ladder 
search  region  whose  track  length  is  equal  to  the  search  track 
length  required  to  complete  an  ideal  ladder  search,  that  is,  a 
track  length  1  =  m«b.   The  degraded  ladder  search  model 
describes  the  result  of  omissions  and  overlaps  in  a  ladder  search 
to  be  such  that  the  probability  of  detection  for  this  random 
search  is  equal  to  the  probability  of  detection  for  the  degraded 
ladder  search.   Consequently,  since  the  area  of  the  ladder  search 
region  is  m»s»b,   for  the  degraded  ladder  search  model: 

(22)  p  =  1  -  exp(-W/s) . 
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Here,  the  requirement  that  the  coverage  factor  W/s  £  1   for 
Equation  21  can  be  relaxed.   However,  it  should  still  be 
considered  as  an  approximate  condition. 

The  condition  that  the  target  be  fixed  within  the 
rectangular  search  region  is  critical  to  both  Equation  21  and 
Equation  22.   However,  these  results  are  also  applicable  to  a 
search  for  a  moving  target  under  the  conditions  that  are 
described  next. 

A  Barrier  Search  Model:   A  target  moves  with  a  constant 
course  and  a  constant  speed  u.   Both  the  target's  course  and  the 
target's  speed  are  known  by  a  searcher.   The  searcher  establishes 
a  barrier  of  width  b  that  is  perpendicular  to  the  target's 
track  and  moves  on  the  barrier  with  a  speed  v  >  u.   The  barrier 
is  designed  so  that  in  a  reference  system  relative  to  the  target 
the  barrier  search  is  a  ladder  search  that  satisfies  the 
conditions  for  a  ladder  search  that  are  given  above.   There  are 
two  cases  to  consider:   1.  The  barrier  is  established  in  front  of 
the  target.   2.  The  barrier  is  established  behind  the  target. 

From  the  search  geometry  for  a  barrier  established  in  front 
of  the  target,  it  can  be  seen  from  Figure  7  below  that 
6  =  sin-1(u/v)  and  d  =  vt  where  r  =  s/(v  +  u)  is  the  time  to 
move  from  one  search  leg  to  the  next.  The  angle  e  and  the 
perpendicular  distance  d  that  depend  on  u,  v  and  s,  and  the 
width  of  the  barrier  b  are  the  quantities  that  are  required  in 
order  to  establish  the  barrier  operationally. 
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Figure  7.  A  barrier  search  track  shown  for  a  barrier  established 
in  front  of  the  target.  The  track  is  shown  in  a  reference  system 
fixed  relative  to  the  earth. 

For  a  barrier  that  is  established  in  front  of  a  target,  one 
of  three  barrier  types  will  result.   A  barrier's  type  is 
determined  by  the  relation  of  the  distance  d  to  the  distance 
g  =  ut  where  the  time  t  =  b/ (v*  -  u2 ) *   is  the  time  to  complete 
a  search  leg  (cross  the  barrier) .   The  barrier  type  is  determined 
as  follows:   1.  For  g  <  d,  the  barrier  is  an  advancing  barrier. 

2.  For  g  =  d,  the  barrier  is  a  stationary  barrier. 

3.  For  g  >  d,  the  barrier  is  a  retreating  barrier. 

For  a  barrier  established  behind  the  target,  there  is  only 
one  barrier  type  and  it  is  called  an  overtaking  barrier.   For  an 
overtaking  barrier,   e  =  sin_1(u/v)   as  for  a  barrier  established 
in  front  of  the  target.   But,  for  an  overtaking  barrier, 
t  =  s/(v  -  u)   and  d  =  vs/(v  -  u)  . 
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Given  the  target  crosses  the  barrier,  the  probability  of 
detection  for  an  ideal  barrier  search  is  given  by  Equation  21 
and  the  probability  for  a  degraded  barrier  search  is  given  by 
Equation  2  2  where  the  terminology  refers  to  the  nature  of  the 
ladder  search  in  the  reference  system  moving  with  the  target. 
A  discussion  of  an  application  of  these  two  equations  to  a  search 
for  a  magnetic  anomaly  target  is  given  in  Reference  15. 
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XII.   A  Target  State  Estimation  Procedure 

A  target  state  estimation  procedure  based  on  bearing 
observations  is  developed  in  this  section  that  generates  point 
estimates  of  a  target's  position  and  velocity  vector  coordinates 
in  a  rectangular  coordinate  system.   The  procedure  is  based  on  a 
model  in  which  bearing  errors  are  unknown  and  are  not  determined 
by  random  variables  with  known  distributions.   Because  of  this, 
confidence  regions  for  the  estimates  are  not  generated  by  the 
procedure.   However,  for  a  moving  target,  it  illustrates  general 
characteristics  of  bearings  only  target  motion  analysis  (TMA) . 
The  model  is  defined  as  follows:   1.  The  target  moves  in  a  plane 
with  a  constant  but  unknown  course  and  speed.   2.  Observations  of 
the  target  are  made  from  known  positions  at  known  times.   3.  The 
observations  provide  only  target  bearings  with  unknown  errors. 
The  model  geometry  is  shown  in  Figure  8. 
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di  =  r-j/sin  (p^   -  Q^) 

[xt(i)  ,YtU)]   estimate 
ri  range  estimate 
j8i  bearing  estimate 


[x0(i) ,y0(i) ]   observer 


x  East 


Figure  8.   The  geometry  of  the  target  motion  analysis  model 
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The  procedure  criterion  is:   For  observations  from  n  positions, 
chose  target  position  estimate  and  target  velocity  component 
estimates  ux  and  Uy  that  make  the  sum  S   of  the  squares  of 
the  algebraic  distance  between  the  estimated  position  and  the 
bearing  line  corresponding  to  the  observed  bearing  a  minimum. 
From  Figure  8,  the  algebraic  distance  can  be  written  as 
di  =  [x(i)  -  x0(i)]-cos  0i   -  [yt(i)  "  y0(i)]*sin  0^.      Because  of 
the  requirement  that  the  target  move  with  constant  course  and 
speed  during  the  encounter,  the  number  of  independent  estimates 
is  reduced  from  2n  to  4:   x^(l),   y^(l),   ux  and  uy.   These 
four  estimates  define  a  target's  motion  by  the  relations: 
xt(i)  =  xt(l)  +  ux.(ti  -  tx)   and  yt(i)  =  yt(l)  +  uy- (ti  -  tx)  . 
To  determine  estimates  of  the  target  state  parameters,  take  the 
partial  derivative  of  the  sum  S  =  E  (d±)2      with  respect  to  each 
of  them.   Then  set  the  four  partial  derivatives  equal  to  zero. 
This  creates  four  linear  equations  in  x^(l)  ,  y-tt1)/  ux  anc*  uy 
whose  solution  are  the  desired  estimates  x^-(l),  y-t(l)/  ux  and 
uv.   In  matrix  notation,  the  equations  can  be  represented  by: 
AX  =  B  where  the  elements  of  X  are  x^(l) ,      yt(l) /   ux»   and 
uy#   A  necessary  condition  for  a  unique  solution  for  X  is  that 
n  >  4.   Otherwise,  the  determinant  of  A  will  be  equal  to  zero. 
The  procedure  can  also  be  used  if  a  target's  course  and  speed  are 
constant  and  known  and,  in  particular,  if  the  target  is 
stationary  so  that  ux  and  uy  both  equal  zero.   In  this  case, 
since  the  number  of  unknowns  is  two,  the  number  of  linear 
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equations  is  also  two  and  a  necessary  condition  for  a  unique 
solution  is  n  >  2. 

Now,  suppose  the  observations  are  at  positions  and  times 
that  correspond  to  the  positions  and  times  of  an  observer  moving 
on  some  constant  course  at  some  constant  speed  (including  zero 
speed) .   In  this  case,  the  observation  position  coordinates  are 
related  by  the  following  equations:   x0(i)  =  x0(l)  +  vx(t^  -  t±) 
and  y0(i)  =  yQ(l)  +  vy(fci  "  fcl)   where  vx  and  vy  are  the 
required  velocity  components  of  the  observer.   Using  these 
equations  of  motion,  the  matrix  equation  AX  =  B  can  be 
transformed  to  the  matrix  equation  AX  =  0  where  the  elements 
of  the  matrix  X  are  related  to  the  elements  of  the  matrix  X 
by  the  equations:   x^  =  xt(l)  -  x0(l)  ,   x2i  =  ytU)  -  Yoi1)  > 

*31  =  ux  "  vx  and  *41  =  uy  "  vy 

Since  the  linear  equations  represented  by  AX  =0  are 

homogenous,  they  do  not  have  unique  solutions  and  consequently 

neither  do  the  equations  represented  by  AX  =  B.   However,  if 

there  is  at  least  one  observation  whose  time  and  position  is  not 

determined  by  the  above  equations  of  motion,  then  the 

transformation  from  X  to  X   cannot  be  made,  and  in  general  a 

unique  solution  for  X  can  be  found.   If  the  observations  are 

made  from  a  platform  that  is  moving  with  a  constant  course  and 

speed,  this  condition  can  be  achieved  by  either  changing  the 

course,  the  speed  or  both  prior  to  completing  the  observations. 

Estimation  models  that  describe  bearing  error  as  a  random 

variable  provide  a  basis  for  determining  confidence  regions  for 
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point  estimates.   A  model  is  developed  in  Reference  16  that  does 
this  for  either  target  bearing  observations  made  from  two  or  more 
points  simultaneously  or  for  a  target  that  is  stationary  relative 
to  the  observation  points. 
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XIII.   Position  Distributions  That  Change  with  Motion 

In  this  section,  target  motion  models  are  discussed 
that  provide  a  basis  for  determining  position  distributions  that 
change  with  time.   In  the  first  model,  a  target  moves  in  a  plane 
with  a  constant  course  and  speed.   In  the  second  model,  a  target 
moves  in  a  plane  with  a  constant  course  and  speed  during  a  time 
step  and  at  the  end  of  the  time  step  changes  either  its  course  or 
speed  or  both. 

Position  Distributions  based  on  the  First  Motion  Model:   Two 
cases  of  this  model  are  considered  that  provide  the  basis  for 
determining  a  position  distribution  analytically.   In  both  cases, 
the  target's  position  coordinates  in  a  rectangular  coordinate 
system  are  determined  by  a  circular  normal  distribution  at  time 
zero.   And,  in  both  cases,  during  the  motion  the  target's 
velocity  does  not  change. 

The  First  Case:   In  this  case,  the  targets 's  coordinates 
X(t)   and  Y(t)   at  time  t  are  determined  by  the  vector  sum 
with  components:   X(t)  =  X(0)  +  Ux*t  and  Y(t)  =  Y(0)  +  Uy»t 
where  X(0)   and  Y(0)   are  the  target's  coordinates  at  time  zero 
and  Ux  and  Uv  are  the  target's  velocity  components.   The 
coordinates  X(0)   and  Y(0)   are  determined  by  a  circular  normal 
distribution  whose  mean  vector  is   (0,0)   and  whose  standard 
deviation  is  a.   And  the  velocity  components  Ux  and  Uv  are 
determined  by  a  circular  normal  distribution  whose  mean  vector  is 
(ttx,«y)   and  whose  standard  deviation  is  au.   The  random  vectors 
[X(0),Y(0)]   and   (Ux»t,Uy»t)   are  independent  and  both  are 
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normal  random  vectors.   Since  they  are  normal  random  vectors,  the 
random  vector   [X(t),Y(t)]   is  also  a  normal  random  vector.   And, 
since  they  are  independent,  its  characteristic  function  is  equal 
to  the  product  of  their  characteristic  functions.   This  implies 
that   [X(t),Y(t)]   has  a  circular  normal  distribution  whose  mean 
vector  is   (ttx«t,ttv«t)   and  whose  variance  is   (a2  +  a^'t2  )  .   For 
a  more  complete  discussion  of  the  basis  for  the  above  argument, 
see  Reference  7. 

The  Second  Case:   The  procedure  that  is  used  to  determine 
the  distribution  of  the  random  vector   [X(t),Y(t)]   for  the 
second  case  can  also  be  used  to  determine  it  for  the  first  case. 
It  is  based  on  the  following  consideration:   In  general,  the 
joint  density  function  of  a  target's  coordinates  at  some  time  t 

greater  than  a  zero  of  time  is  determined  as  follows: 

Ooo 
-oofX(0)  ,Y(0)  (q»s'°)  fv,w(v,w)  dvdw 
where  q  =  x  -  vt  and  s  =  y  -  wt  and  V  =  Ux,  W  =  Uv,  v  =  ux 
and  w  =  Uy.   Equation  23  can  be  developed  as  follows:   To  first 
order,  fx(t)  Y(t) (x»Y't) • Sx'Sy     is  the  probability  that  a 
target's  coordinates  are  in  an  element  of  area  6x*6y     and  for 
given  values  of  v,  w  and  t,  target  positions  in  an  element  of 
area   <Sq«<5s  will  be  translated  to  an  element  of  area  Sx*Sy 
that  is  identical  in  form  and  size  to  6q«<Ss.   And  so,  to  first 

order,  fx(0)  *  Y(0)  (3/s'°)  *  Sx'  SV  fV,w(u' v)  *  5v* <Sw   ^s  tne 
probability  that  the  target's  coordinates  at  time  0  are  in  an 
element  of  area   6q»<Ss  that  is  located  such  that  the  target's 
coordinates  will  be  in  the  element  of  area   <5x»6y  at  time  t 
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since  x  =  q  +  vt  and  y  -  s  +  wt.   And,  to  first  order,  the 
sum  of  such  probabilities  for  all  pairs  of  values  of  v  and  w 
is  also  the  probability  that  the  target's  coordinates  at  time  t 
are  in  the  element  of  area  Sx'Sy.      In  the  limit  after  equating 
the  two  expressions  for  this  probability  and  cancelling  the 
common  factor  6x»6y,    Equation  23   results. 

In  the  first  case,  the  joint  distribution  of  the  velocity 
components   (UX,UV)   is  circular  normal  with  mean  vector 
(ux,Uy)   and  standard  deviation  au  and  the  joint  distribution 
of   [X(0),Y(0)]   is  circular  normal  with  mean  vector   (0,0)   and 
standard  deviation  a.      In  the  second  case,  the  target's  speed 
u   is  known  but  its  course  is  equally  likely  to  have  any  value 
between  0  and  2ir ,  so  Ux  =  u»sin  *  and  Uv  =  u»cos  *  where  * 
is  a  random  variable  with  a  uniform  distribution  and  values 
0  <  <p  <   2w .      In  this  case,  only  the  single  random  variable  *  is 
required  to  determine  the  target's  velocity.   With  the  velocity 
components  in  polar  coordinates,  the  integral  of  Equation  23   is 
a  single  integral  over  <p     and  the  integrand  of  the  integral  is 
(l/2iro2  )  »exp[-(q2  +  s2  )/2a2  ]  •  (1/2tt)    where  now  q  =  x  -  u»t»sin  <p 
and  s  =  y  -  u»t*cos  <p.      The  result  of  the  integration  is: 
(24)  (l/27ra2)  exp{-[x2  +  y2  +  (u»t) 2  ]/2a2  }  I0[(x2  +  Y2  )  *  *u»t/a2  ] 
where  I0   indicates  the  hyperbolic  Bessel  function  of  zeroth 
order.   Expression  24   is  plotted  in  Reference  6   for  several 
values  of  t.   There,  a  single  variable  r  is  substituted  for 
(x2  +  y2)"5,  the  the  target's  range  from  the  origin.   The  plot  of 


52 


the  expression  shows  a  limiting  characteristic  of  the 
distribution  that  can  be  indicated  as  follows:   First  replace 
(x2  +  y2 ) i   by  r   in  Expression  24.   Then,  by  first  multiplying 
and  then  dividing  the  resulting  expression  by  exp(-r»u»t/a2 ) ' 
the  following  expression  can  be  obtained: 

(25)  l/(2na2  )exp{-[l/ (2a2  )  ]  (r  -  u»t) 2  }I0  (r»u-  t/cr2  )  exp(-r»u»t/a2  )  . 
The  second  factor  in  Expression  25   is  proportional  to  the 
density  function  of  a  normal  random  variable  whose  mean  is  u«t 
and  whose  standard  deviation  is  a.      And  the  last  two  factors  are 
such  that  their  product  changes  slowly  as  a  function  of  r 
relative  to  the  second  factor.   The  consequence  of  this  is  that  a 
plot  of  Expression  25   against  r   for  values  of  t  greater  than 
4»a/u  has  the  appearance  of  a  normal  density  function.   In 
particular,  this  indicates  that  u»t   in  the  exponent  of  the 
second  factor  in  Expression  25  can  be  viewed  as  the  radius  of  an 
average  furthest  on  circle. 

A  target's  rectangular  coordinates  and  its  range  and  bearing 
from  the  origin  are  related  by:   X  =  R  sin  6  and  Y  =  R  cos  8. 
By  using  these  equations,  Expression  24   can  be  transformed  to: 

(26)  (l/2?r)  (r/27ra2  )exp{-[r2  +  (ut) 2  ]/2a2  }I0  (rut/a2  ) 

which  is  the  joint  density  function  of  the  random  variables  R 
and  9.   The  marginal  density  function  of  R  can  be  obtained 
from  Expression  26  by  integrating  this  joint  density  function 
over  the  possible  values  of  0  which  in  radians  is  over  the 
interval   0   to  2w.   Consequently,  the  marginal  density  for  R 
can  be  obtained  from  Expression  26  by  multiplying  it  by  2tt. 
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Tabulated  values  of  the  cumulative  distribution  function  FR(r;t) 
for  the  marginal  distribution  of  R  for  the  second  case  are 
listed  in  Reference  17. 

Position  Distributions  Based  on  the  Second  Motion  Model:   In 
the  second  motion  model  a  target's  motion  is  not  restricted  to  a 
constant  course  and  speed  throughout  the  motion.   At  the  end  of  a 
time  step,  a  target's  course  or  speed  or  both  can  change.   In 
general,  this  leads  to  a  position  distribution  that  can  not  be 
described  analytically.   And,  in  general,  a  monte  carlo 
simulation  method  is  required  in  order  to  determine  a  position 
distribution  that  is  based  on  such  a  model.   As  an  example  of  a 
case  in  which  the  distribution  can  be  described  analytically,  see 
Reference  18.   As  an  example  of  a  case  in  which  it  can  not, 
suppose  a  target's  initial  position  is  described  in  terms  of  a 
number  assigned  to  a  subregion  in  the  xy-plane  where  the  number 
assigned  represents  the  probability  that  the  subregion  contains 
the  target  at  an  initial  time.   In  addition,  suppose  for  each 
subregion  a  course  and  speed  distribution  is  determined  by 
assigning  numbers  to  course  and  speed  pairs  where  a  number 
represents  the  probability  the  target  will  have  the  course  and 
speed  at  the  initial  time  given  it  is  in  the  subregion  at  that 
time.   Next  suppose  for  each  course  and  speed  pair  there  is  a 
time  distribution  that  determines  the  duration  of  the  course  and 
speed  pair  and  that  the  time  distribution  is  determined  by  a 
number  assigned  to  each  discrete  time  point  where  the  number 
represents  the  probability  that  the  target's  course  and  speed 
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pair  will  be  determined  by  a  new  course  and  speed  distribution. 
By  extending  this  kind  of  procedure  and  then  implementing  it  in  a 
monte  carlo  simulation,  one  can  generate  complex  position 
distributions  that  describe  a  target's  position  at  discrete  time 
points  in  terms  of  probabilities  assigned  to  subregions  of  a 
region  in  which  the  target  is  contained. 
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XIV.   Position  Distributions  That  Change  with  Search 

Suppose  a  target's  position  at  some  time  is  described 
by  a  position  distribution.   Now  suppose  information  becomes 
available  that  a  search  has  been  conducted  for  the  target  and 
that  the  target  has  not  been  detected.   Or  suppose  the 
information  is  that  the  target  has  been  detected.   In  the  first 
case,  negative  information  is  available  that  can  be  used  to 
modify  the  position  distribution.   In  the  second  case,  positive 
information  is  available  that  can  be  used  to  modify  the  position 
distribution.   In  both  of  the  cases,  a  target's  position 
distribution  is  assumed  to  be  specified  in  terms  of  a  set  of 
discrete  probabilities  where  each  probability  corresponds  to  a 
subregion  of  the  region  that  contains  the  target  and  each  is  the 
probability  that  the  target  is  in  that  subregion. 

Position  Distributions  and  Negative  Information:   For  a 
region  that  contains  a  target  and  consists  of  n  subregions,  let 
the  event  S^  =  {the  target  is  in  the  ith  subregion}.   And  let 
the  event  C  =  {no  contact).   Then,  given  no  contact  in  a  search 
of  the  region,  the  targets' s  position  distribution  can  be 
modified  as  follows: 

(27)    P(Si|C)   =   PfClSiJ.PfSiJ/PfC) 

where  i  =  1,2,  •••  ,n  and  P(C)  =  Z  P(c| Sj ) »P(Sj )   with  the  sum 
from  j  =  1,2,  •••  ,n.   Note  that  Equation  27   can  be  obtained  by 
using  Bayes  theorem.   To  illustrate  how  Equation  26  might  be 
used,  suppose  that  a  search  in  a  subregion  is  considered  to  be  a 
random  search  and  that  the  sweep  width  of  detection  system 
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against  the  target  depends  on  subregion  being  searched.   For  this 
case,  let  A^  be  the  area  of  the  ith  subregion  and  let  W^  be 
the  sweep  width  in  that  subregion.   Then,  given  no  contact  in  a 
search  of  a  subregion,   P(c|Si)  =  exp[-(Wj/  ljJ/AjJ   where   1^   is 
the  track  length  of  the  searcher  in  the  itn  subregion.   Given 
values  for  P(SjJ  ,  6±     and  1^  for  i  =  1,2,  •••  ,n  a  position 
distribution  can  be  determined  that  has  been  modified  by  the 
negative  information. 

Position  Distributions  and  Positive  Information:   In  the 
case  of  positive  information,  the  event  C  =  {a  contact}   occurs. 
This  event  is  the  union  of  two  mutually  exclusive  and  exhaustive 
events:   Tc  =  {a  true  contact)   and  Fc  =  {a  false  contact). 
Relative  to  the  Venn  diagram  of  Figure  1,  the  event  S^  defined 
above  corresponds  to  H^,   C  corresponds  to   D^,   Tc 
corresponds  to   (Dj_  n  H^)   and  Fc  corresponds  to   (D^  n  Hq)  . 
For  generality,  suppose  true  contacts  do  not  localize  a  target  to 
a  single  subregion.   Then  after  a  search  of  a  region  that  has 
resulted  in  a  contact,  the  target's  position  distribution  can  be 
modified  as  follows: 

(28)  P(Si|C)  =  P(Si|Tc)-P(Tc|C)   +   P(Si|Fc)-P(Fc|C) 
where  i  =  1,2,  •••  ,n  since  Tc  =  Tc  n  C,   Fc  =  Fc  n  C  and 
P(SjJC)  =  (P[(Si  n  Tc)  +  P(Si  n  FC)])/P(C).   The  probability 

p  =  P(TC|C)   has  been  called  the  credibility  of  the  contact.   In 
terms  of  p  Equation  28  becomes: 

(29)  P(Si|C)  =  P(Si|Tc).p  +  P(Si|Fc)-(l  -  p). 
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For  the  model,   P(SjjTc)   can  be  determined  by  the  coverage 
characteristics  of  the  detection  system  that  is  used  to  make  the 
contact.   In  particular,   P(SjJtc)  =  P(TC| S^) •P(Si)/P(Tc)   where 
P(TC)  =  S  P(Tc|Sj) •P(Sj)   and  the  sum  index  j  =  1,2,  •••  ,n. 
Then,  using  the  correspondence  between  the  events  S^   and  H^ 
and  the  events  Tc  and   (D^  n  H^  ,   a  definition  for  P(Tc|Si) 
is  P(Tc|Si)  =  (Pd) i  where   (Pd) i  is  an  average  probability  of 
detection  over  the  i*-*1  subregion.   In  one  positive  information 
model,   P(SjjF)  =  P(SjJ   (false  contacts  supply  no  information 
about  a  target's  location),  and  p  is  determined  subjectively 
based  on  factors  associated  with  the  detection  system  used  to 
make  the  contact.   The  probability  determined  by   1  -  p  has  been 
called  the  false  alarm  probability.   However,  it  is  not  the 
probability  pf  =  PfD^jHo).   Since   1  -  p  =  P(Fc|c),  it  is  the 
probability  PCHqID!). 

The  above  approach  can  also  be  used  to  develope  a  procedure 
to  modify  a  position  distribution  given  negative  information.   To 
do  this,  note  that  the  event  C  is  the  union  of  two  mutually 
exclusive  and  exhaustive  events:   T^  =  (a  true  no  contact)   and 
Fq   =  {a  false  no  contact).   Relative  to  the  Venn  diagram  of 
Figure  1,   C  corresponds  to  D0,   T5  corresponds  to   (Dq  n  Hq) 
and  F^  corresponds  to   (Dq  n  HJ .   In  analogy  to  Equation  29, 
(30)   P(Si|C)  =  P(Si|Tc)'P  +  P(Si|Fc)-(l  "  P) 

where  p~  =  PCT^Ic)  .   Relative  to  the  Venn  diagram  of  Figure  1, 
p  =  P(D0|H0).   In  analogy  to  the  choice  for  P(SjJtc)  above,  one 
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could  choose  PfSjjTc)  =  P(Si)   and  P(S-jjFc)   could  be 
determined  by  the  coverage  characteristics  of  the  detection 
system.   In  particular,   PfSjjFc)  =  P(Fc  I  S±)  •PfS-jJ/Fc)   where 
P(Fq)  =  S  P(FclSj) »P(Sj)   and  the  sum  index  j  =  1,2,  ••  ,n. 
Using  the  correspondence  between  the  events   F5  and  D0  n  H^ 
one  can  write  P(Fcls^)  =  1  -  (Pa)i. 

Both  in  the  case  of  this  procedure  for  negative  information 
and  the  procedure  for  positive  information,  (Pd) i   can  be  viewed 
as  a  measure  of  the  specified  performance  of  a  detection  system 
and  both  p  and  p   can  be  viewed  as  measures  of  the  degradation 
in  its  specified  performance  do  to  various  operational  factors. 
In  particular,  for  p  =  0,  the  expression  for  P(SjJc)   is  the 
same  for  the  two  negative  information  procedures  if   P(C|Si)  is 
determined  in  the  same  way  for  both  cases.   In  the  view  just 
expressed,  this  implies  that  the  first  procedure  applies  to  a 
search  in  which  a  detection  system  is  operated  so  that  it 
achieves  its  specified  performance.   In  the  case  of  the  positive 
information  procedure,  if  p  =  1,  then  the  expression  for  P(SjJc) 
is  the  same  as  that  which  can  be  obtained  from 
(31)    P(Si|C)  =  P(C|si)«P(Si)/P(C) 

where  P(C)  =  2  P(C| Sj ) *P(Sj )   and  the  sum  index  j  =  1,2,  •••  ,n 
if  P(C|Si)  is  determined  in  the  same  way  for  both  cases. 
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XV.   Search  Models  and  Search  Theory 

Search  theory  provides  a  basis  for  determining  optimal 
search  plans  for  a  target  whose  motion  and  location  are 
determined  within  some  bounds.   Here,  an  optimal  search  plan  is 
one  for  which  the  probability  of  finding  a  target  within  a  given 
length  of  time  is  a  maximum,  the  expected  time  to  find  a  target 
is  a  minimum  given  the  target  is  found  or  a  search  plan  for  which 
some  other  optimal  search  criterion  is  satisfied. 

Search  theory  results  are  based  on  models  of  the  search 
process.   To  the  degree  that  a  search  model  describes  a  search 
process,  an  optimal  search  plan  for  a  target  that  is  based  on  the 
search  model  should  provide  guidance  for  the  development  of  an 
operationally  feasible  search  plan.   However,  because  of  the 
limitations  of  analytical  search  models,  an  optimal  search  plan 
that  is  based  on  an  analytical  search  model  may  give  only  initial 
guidance  in  this  regard.   The  optimal  search  plans  that  are 
described  below  illustrate  this.   The  search  plans  are  based  on 
the  random  search  model.   Because  of  this,  the  requirement  on  the 
location  of  search  track  segments  is  not  realizable  and  the  time 
to  resolve  false  alarms  is  ignored. 

Optimal  search  plans  based  on  search  models  implemented 
through  a  monte  carlo  simulation  are  not  considered  here. 
However,  with  sufficient  information,  such  plans  have  the 
potential  of  being  both  implementable  and  more  optimal  in  a  real 
sense  than  an  optimal  search  plan  based  on  an  analytical  search 
model . 
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Three  Optimal  Search  Plans:   The  three  optimal  search  plans 
differ  through  their  definition  of  optimality.   However,  each  one 
is  based  on  the  following  search  model:   A  target  is  fixed  at 
some  point  in  a  region  that  consists  of  n  subregions.   A  search 
in  a  subregion  is  a  random  search  in  the  sense  of  the  definition 
in  Section  X  and  a  searchers  sweep  width  there  is  a  constant. 
In  addition,  a  search  of  a  subregion  will  not  detect  a  target 
which  is  in  another  subregion.   To  determine  a  plan,   let 
S^  =  {the  target  is  in  subregion  i}   for   i  =  1,2,  •••  ,n  and 
let  Pi  =  P(Si)   be  the  prior  probability  that  the  target  is  in 
the  itn  subregion.   Let  W^  be  the  sweep  width  in  the  itn 
subregion.   Let  6^  -   Aj/W^  where  A^   is  the  area  of  the  itn 
subregion  and   6^   is  the  expected  track  length  to  find  the 
target  by  a  search  of  the  itn  subregion  given  the  target  is  in 
the  itn  subregion,  a  characteristic  length.   The  probability  P 
that  the  target  will  be  detected  by  a  random  search  is  given  by: 
(32)      P  =  S  [1  -  exp(-li/6i)].pi 

where  the  sum  index   i  =  1,2,  •••  ,n   and   1^   is  the  track 
length  of  the  search  in  the  itn  subregion. 

The  first  optimality  criterion  is:   Choose   1^   so  that  P 
is  a  maximum  subject  to  the  two  constraints:   1.  1  =  2  lj_   and 
2.  li  >  0  where  the  index   i  =  1,2,  •••  ,n.   Determining  this 
choice  is  a  nonlinear  optimization  problem  whose  solution  is 
given  in  Reference  19.   It  is: 

l*/6i  =  ln(pi/6i)  -  L(k)    i  =  1,2,  •••  ,  k 

(33) 

li/Si   =0  i  =  k+l,k+2,  •••  ,n 

61 


where  L(k)  =  (1/Z6j )  *Z[6j  «ln(pj/6j )  ]  +  l/S<Sj   and  the  sum  index 
j  =  1,2,  •••  ,  k,  where  the  subregions  are  relabeled  so  that  the 
following  order  relation  holds:   Pi/<*i  >  P2/*2  >  •••  >  Pn/*n  and 
where  k  is  chosen  so  that  for  k+1  the  solution  for  lfc+i 
using  L(k+1)   is  either  negative  or  zero. 

The  second  optimal ity  criterion  is:   Choose   1^   so  that  P 
is  a  maximum  subject  to  the  two  constraints:   1.  c  =  2  c^  and 
2.  c^  >  0  where  the  index  i  =  1,2,  •••  ,n,  c-^  =  fcj/li   is  the 
cost  of  the  search  in  the  ith  subregion  and  k^   is  the  cost  per 
unit  track  length  in  that  subregion.    For  this  criterion,  the 
solution  to  the  corresponding  nonlinear  optimization  problem  can 
be  obtained  from  Equation  33   by  replacing   6^  by   €^  =  kj/$i 
and  labeling  the  subregions  so  that  Pi/^i  >  P2/€i  >  •••  >  Pn/€i* 
The  basis  for  this  can  be  seen  by  replacing  li/<$i  by   its 
equivalent  ci/£±      in  the  exponential  term  in  Equation  32. 

The  third  optimal ity  criterion  is:   Choose   1^   so  that  the 
expected  utility  of  finding  the  target  is  a  maximum  subject  to 
the  two  constraints:   1.   1  ■  E  lj  and  2.   1^  >  0  where  the 
index  i  =  1,2,  •••  ,n.   For  this  criterion,  the  solution  to  the 
corresponding  nonlinear  optimization  problem  can  be  obtained  from 
Equation  31  by  replacing  p^  by  q^  where  q^  =  ui*Pi  and  ui 
is  the  utility  of  finding  the  target  given  it  is  in  the  itn 
subregion.   And,  in  addition,  labeling  the  subregions  so  that 
3l/5l  >  <32/<^2  >  •  •  •  >  <5n/^n*   Tne  basis  for  this  can  be  seen  by 
multiplying  the  summation  term  in  Equation  32  by  u^  so  that 
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the  resulting  equation  gives  the  expected  utility  of  finding  the 
target  with  the  search. 

Equation  33   can  be  used  to  determine  an  order  of  search  for 
the  subregions  which  will  effectively  minimize  the  expected  track 
length  required  to  detect  a  target  given  it  is  detected.   To  do 
this,  divide  the  available  track  length  1   into  units  small 
enough  so  that  with  a  single  unit  only  the  1st  subregion  would  be 
searched.  Then  allocate  one  unit  to  the  search  of  the  1st 
subregion.   If  the  search  is  unsuccessful,  determine  the  optimum 
allocation  for  two  units.   Then  search  with  a  second  unit  so  that 
the  first  search  with  the  first  unit  plus  the  second  search  with 
the  second  unit  satisfy  the  optimum  allocation  for  two  units.   If 
the  search  is  unsuccessful,  continue  in  this  fashion  until  either 
the  target  is  found  or  all  the  track  length  is  expended.   That 
this  allocation  order  will  effectively  minimize  the  expected 
track  length  required  to  detect  a  target  given  it  is  detected  can 
be  argued  as  follows:   Let  L  be  the  track  length  at  detection, 
let   lu  be  a  unit  of  track  length  and  let  n  be  the  number  of 
units.   Then  the  value  of  the  probability  P(L  <  i'lu)   that  the 
target  will  be  detected  on  or  before  the  i*-*1  step  of  the  search 
for  the  given  allocation  order  will  be  greater  than  or  equal  to 
its  value  for  any  other  allocation  order  with  the  same  allocation 
step  size.   Since  the  value  of  P(L  <  1)   will  be  equal  to  its 
value  for  any  other  allocation  order  of  the  optimum  allocation 
and  since   P(L  <  i«lu|L  <  1)  =  P(L  <  i-lu)/P(L  <  1),  the  value  of 
the  distribution  function  FL(i»lu|L  <  1)  =  P(L  <  i»lu|L  <  1) 
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will  be  greater  than  or  equal  to  its  value  for  any  other 
allocation  order.   This  implies  that  the  expected  track  length 
given  detection  E(l|l  <  1)  =  S  [1  -  FL(i«lu|L  <  1)]   where  the 
sum  index  1=1,2,  •••  ,n  is  effectively  a  minimum  for  the 
given  allocation  order.   A  search  based  on  the  optimum  allocation 
given  by  Equation  33   and  the  given  allocation  order  is 
equivalent  to  the  following  search:   After  an  allocation  of  track 
length   lu  and  an  unsuccessful  search,  new  values  for  P(S^)   are 
calculated  using  Equation  27   and  then  Equation  3  3   is  used  with 
these  new  values  to  determine  the  next  optimum  allocation.   A 
discussion  of  this  procedure  is  given  in  Reference  6.   And  an 
example  of  its  application  is  given  in  Reference  20. 

Equation  33   also  defines  an  optimal  search  plan  for  a 
detection  system  that  searches  beams  and  can  be  described  by 
Equation  31  by  replacing   1^  by  t^  where  t^   is  the  time  the 
i"1  beam  is  searched  and  by  replacing  6±     by  r^  where  r^,   a 
characteristic  time,  is  the  expected  time  to  detect  the  target  by 
a  search  of  the  itn  beam  given  the  target  is  in  the  itn  beam. 

For  a  more  extensive  discussion  of  search  theory  and  its 
application  to  military  operations  research,  see  Reference  21. 
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